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Abstract 

We consider a very general dilaton-axion system coupled to Einstein-Hilbert gravity in ar- 
bitrary dimension and we carry out holographic renormalization for any dimension up to and 
including five dimensions. This is achieved by developing a new systematic algorithm for itera- 
tively solving the radial Hamilton- Jacobi equation in a derivative expansion. The boundary term 
derived is valid not only for asymptotically AdS backgrounds, but also for more general asymp- 
totics, including non-conformal branes and Improved Holographic QCD. In the second half of 
the paper, we apply the general result to Improved Holographic QCD with arbitrary dilaton 
potential. In particular, we derive the generalized Fefferman- Graham asymptotic expansions 
and provide a proof of the holographic Ward identities. 
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1 Introduction 

Understanding the holographic dictionary for holographic models in non asymptotically AdS spaces 
has been a long standing problem. It has been a pressing question ever since physically promis- 
ing holographic dualities involving non asymptotically AdS backgrounds, such as the Klebanov- 
Strassler PQ and Maldacena-Nuhez [2j backgrounds, were found, but it has become even more 
relevant with the recent interest in the phenomenological application of holography to condensed 
matter physics and models of QCD. Even though numerous attempts have been made to understand 
aspects of the dictionary of some of these systems, there are very few cases where a systematic and 
extensive understanding has been achieved for non asymptotically AdS backgrounds. These include 
the analyses of the dictionary for non-conformal branes (3[ [4] and for Schrodinger backgrounds [5] . 

However, it cannot be overemphasized that the process of understanding the holographic dictio- 
nary, in any holographic model and any background, even beyond the supergravity approximation, 
can be split into two conceptually distinct steps. The first step is intrinsically related with the 
"bulk" holographic model. Namely, one must identify a suitable boundary in the bulk theory and 
construct a reduced phase space of the theory in terms of data on that boundary. This step is 
exactly analogous to the Fefferman-Graham reconstruction of the bulk geometry in asymptotically 
hyperbolic manifolds from boundary data [6 J . A systematic way of addressing this question and its 
connection to a certain variational problem at infinity in the most general setting was discussed in 
[7]. The approach developed in [7] in principle allows one to algorithmically construct this reduced 
phase space for any bulk model. Having completed this step, one not only has achieved a reformula- 
tion of the bulk dynamics in terms of a symplectic space of boundary data that can unambiguously 
be identified with the symplectic space of renormalized observables in any holographically dual 
theory, but also has automatically made the variational problem of the bulk theory well defined, 
which implies that the on-shell action is finite [HI \7\ . 

Only once this first step has been completed, one can directly compare the symplectic space 
of boundary data with the symplectic space of gauge-invariant observables in any candidate holo- 
graphic dual. This mapping is simply the classical version of the Hilbert space isomorphism one 
expects in a fully quantum mechanical holographic duality. On the bulk side one has the symplectic 
space of a classical system, being that classical strings or classical gravity, while on the field theory 
side the Hilbert space reduces to a classical symplectic space in some limit where the number of 
degrees of freedom becomes infinite [9]. 
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In this paper we consider a generic dilaton-axion system coupled to Einstein-Hilbert gravity in 
arbitrary dimension with the action (|2.ip . This system contains the standard dilaton-axion system 
in AdSs dual to the complexified coupling of J\f = 4 super Yang-Mills in four dimensions, non- 
conformal branes OH], as well as Improved Holographic QCD [10] as special cases. The last two 
examples admit non asymptotically AdS vacua, and so the standard dictionary for asymptotically 
AdS gravity is not applicable. Our aim here will be to carry out this two-step procedure outlined 
above for this general dilaton-axion system and to explore in more detail the consequences for the 
model of Improved Holographic QCD. 

The semi-phenomenological holographic model dual to large iV c Yang-Mills theory put forward 
in [lOj (see [TTJ for an extensive review) is based on the five dimensional two-derivative (Euclidean) 
bosonic supergravity action 

S = -M^N 2 , J Sx^fg (R[g] - fX^d^X - Z(X)d flX d fl X + V(X)) , (1.1) 

where the Planck mass = l/g^l^ is related to the five dimensional Newton's constant by 
(167rG5) _1 = M^N^. The field content of this action consists of the five dimensional metric g^ u , 
a dilaton A, and an axion \- These are respectively designed to describe the dynamics of the 
lowest dimension gauge-invariant operators of pure Yang-Mills theory, namely the stress tensor, 
Tij, Tr(F 2 ) and Tr(.F'-F). In particular, A is proportional to the 't Hooft coupling, N c g YM , while 
X is related to the instanton angle 9ym- The constants of proportionality are not known a priori, 
but can be determined by comparing the perturbative UV expansion of the beta functions for the 
't Hooft coupling and 9ym with the corresponding holographic beta functions for the bulk fields 
A and x respectively. These proportionality constants are nevertheless scheme dependent and so 
they do not affect the value of any physical observable [11] . 

The holographic model is defined by the potential V(X) and the function Z(A), as well as the 
constant £ 7^ 0, corresponding to the normalization of the kinetic term of the dilaton|^| Anticipating 
that such a model could possibly originate in a non-critical string theory in five dimensions [10J, 
the metric and the dilaton are expected to come from the NSNS sector while the axion comes from 
the RR sector. As was argued in [10], this implies that the kinetic term of the axion should be 
0(1/N^) relative to the Einstein-Hilbert term, the dilaton kinetic term, and the scalar potential. 
Hence, Z(X) = C(l/iV c 2 ), while f and V(X) are <D(N°). 

The form of the functions Z(X) and V(X) can be constrained by physical input. Firstly, asymp- 
totic freedom means that at the UV the theory is conformal and so the dual string vacuum should 
be asymptotically AdSs. Of course, it also means that the theory is weakly coupled at the far 
UV, and so a two-derivative gravity approximation of the form (jl.ip cannot a priori be trusted 
in that limit. However, the assumption in [TO] is that the effect of the higher derivative terms 
affecting the dynamics of the lowest lying fields can somehow be incorporated into a an effective 
cosmological constant in the potential U(A), thus allowing the two-derivative action (11. lj) to admit 
asymptotically AdS solutions. Although this argument is admittedly difficult, if not impossible, 
to defend from a string theoretic point of view, the attitude in [10J is that this assumption can 
be a posteriori justified by the success of the model in reasonably describing various qualitative 
properties of Yang-Mills theory. Accepting this assumption, means that the dilaton A must vanish 
in the far UV and the functions V(X) and Z(X) must admit Taylor expansions of the forrrJl 

/ 00 \ 00 

v w = p 1 + E v - xn . z w = (K^r 1 E (i-^) 

V n.= l / n=0 

where i is the radius of the AdS corresponding to the UV fixed point, and V n and Z n are 0(N®). The 
original motivation for these expansions was that they should be the holographic analogue of the 

1 In [10] £ 2 = 4/3, but here we prefer to keep it arbitrary. 

2 Potentials that contain non-analytic terms at order higher than C9(A 2 ) have been considered in the literature. 
The analysis below applies to such more general potentials as well. 
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perturbative expansions of the beta functions f3\ and f3 x of the 't Hooft coupling and instanton angle 
respectively [10]. In particular, the coefficients V% and V% were argued to be related respectively to 
the one- and two-loop beta function of the 't Hooft coupling. In order to account for the logarithmic 
running of the coupling, therefore, one demands that V\ 7^ 0. As we shall see below, the requirement 
that V\ 7^ is indeed crucial for the asymptotic solutions to logarithmically deviate from strictly 
AdS asymptotics. However, the beta function is a scheme-dependent quantity which does not 
correspond to any physical observable. As we shall see below, the holographic Ward identities only 
involve physical, renormalization group invariant quantities. 

Further constraints on the functions V(X) and Z(X) are imposed by the IR properties of the 
model, i.e. as A — > 00. Specifically, confinement and the absence of certain 'bad' singularities in 
the IR require that [10] 



V(X) ~ A 2 « (log A) p , as A - 00, with ^ < ? f^> (1.3) 

Moreover, requiring an asymptotically linear glueball spectrum, m 2 ~ n, uniquely picks out Q = 
2/3, P = 1/2. Certain conditions for the strongly coupled limit of Z(X) are also necessary in some 
cases [101 . Here we will keep the functions V(X) and Z(X) completely general, however, assuming 
only thalU 

V(X) = d{d ~ l) (1 + V\X + V 2 X 2 ) + V(X), V 1 ^0,V(X) = o(X 2 ),asX^0. (1.4) 

In particular, we allow V(X) to contain non analytic terms provided they are subleading compared 
to A 2 as A — > 0. 

Given the functions V(X) and Z(X), the Yang- Mills vacuum is described by extrema of (jl.ip 
with four-dimensional Poincare invariance. Such solutions are domain walls of the form 

ds 2 = dr 2 + e 2A(r) dx i dx i , X = A(r), x = x(r)- (1-5) 

Such backgrounds are extrema of (jl.ip provided they satisfy 

A 2 - -^—-^ (£ 2 A- 2 A 2 + Z(X)x 2 + V(\j) = 0, 
A + dA 2 -j^jV(X)=0, 

2fX- 2 (x + dAX - A^A 2 ) - Z{X) X 2 + V'(X) = 0, 

Z(X) X + Z'(X)Xx + dZ(X)A x = 0, (1.6) 

where ' denotes differentiation w.r.t. the radial coordinate r. These equations are automatically 
solved provided A, X and x satisfy the first order flow equations 

A—jI^x), A = r^^, i-z-HX)?^, (1.7) 

where the 'superpotential' W(X, x) is determined by the scalar potential V(X) and the function 
Z(X) via the equation 



3 We will consider an arbitrary boundary dimension d from now on, in order to maintain a wider applicability of 
our results. 
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In particular, any solution W(\, %) of (jl.8p uniquely specifies a Yang-Mills vacuum. Notice that 
the last equation in (|1.6p can be integrated exactly to obtain 



X = cZ~ l (\)e 



-dA 



(1.9) 



where c is an integration constant. From the last equation in (jl.7p then follows that 



dW(\, X ) 
dx 



= ce 



-dA 



(1.10) 



As we will see later, this result has very significant consequences for the form of the counterterms 
required to make the variational problem for the action (jl.ip well defined. As far as possible 
vacuum solutions are concerned, this relation means that there are two broad classes of vacua 
depending on whether c is zero or not. Namely, if c = 0, then W is independent of the axion and 
so the axion is just a constant in the background, corresponding to the value of the instanton angle 
Qym- Vacua with c ^ allow for a non-trivial axion profile, corresponding to giving a VEV to 
the operator Tr(FF). In particular, using the result (|3.44p below, and assuming O x = Tr(FF), 
we obtain (Tr (FF)} = c/k 2 . Additionally, generic vacua are classified according to whether the 
operator dual to the dilaton acquires a vacuum expectation value. As follows from the exact one- 
point function (|3.40p . this happens iff the UV expansion of the superpotential W(X,x) contains a 
non-zero term of the form e~ dA ~ \ du exp(— d/b^X), where 



Poincare invariant vacua are, therefore, parameterized by three constants: the dilaton and axion 
VEVs, as well as the instanton angle 9ym- As we shall see, the source of the dilaton is a gauge 
freedom and its value can be thought of as the energy scale. In particular, it's value is not a physical 
observable and does not correspond to a coordinate in the moduli space of vacua. This is in good 
agreement with what one expects from a holographic model describing pure Yang-Mills theory. 

The rest of the paper is organized as follows. In Section [2] we define the general dilaton- axion 
gravity theory we will consider and we formulate its dynamics in terms of a radial Hamiltonian. 
We then develop a systematic iterative procedure for solving the Hamilton-Jacobi equation in a 
derivative expansion for this general class of theories. Using this procedure the full boundary term 
that makes the variational problem at infinity well defined is derived for spacetime dimension up 
to and including five dimensions. This result, summarized in Tables [2] and [31 is the main result of 
the paper. In Section [3] we apply this general result to IHQCD. After writing down explicitly the 
boundary counterterms for IHQCD, we systematically derive the generalized Fefferman-Graham 
asymptotic expansions and give explicit expressions for the exact renormalized one-point functions 
in the presence of sources for the stress tensor, dilaton and axion operators. These are given in 
equations (|3.4ip . (|3.40p and (|3.44|) respectively. The section ends with a detailed discussion of the 
asymptotic bulk diffeomorphisms that preserve the form of the asymptotic expansions, leading to a 
proof of the holographic Ward identities. Some concluding remarks follow in Section HI while some 
technical details of the main calculation are presented in Appendix EJ Finally, in Appendix iBl we 
apply the main result to the case of a constant dilaton potential, and derive general expressions for 
the exact one-point functions in the presence of sources for the fully coupled stress tensor-dilaton- 
axion sector of N = 4 super Yang Mills in four dimensions. To our knowledge this is the first time 
that the general form of these one-point functions has been derived. 

2 Boundary term for generic dilaton-axion gravity 

Since the action (|l.ip already contains essentially arbitrary functions of the dilaton, deriving the 
appropriate boundary term that makes its variational problem well defined is not much easier than 



2£ 2 1 (d - 1)V 2 




v = 



d-l d £2fc2 
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considering instead the slightly more general action 

S = ~^{J M dd+lx ^ (RM ~ d^d^ - ZMS^X + V(<p)) + J d d x^2K^ , (2.1) 

where the spacetime dimension is taken to be arbitrary and we have introduced a canonical dilaton 
field, <p, related to the dilaton A in (jl.ip by tp = £ log A. Moreover, we have added the standard 
Gibbons-Hawking term [12] and the constant k is related to Newton's constant in d + 1 dimensions 
by k 2 = MG d+ i. 

Note that this action contains as special cases a very large class of theories considered in the 
literature. Apart, from IHQCD, other special cases include the standard dilaton-axion of N = 4 
super Yang-Mills^], as well as non-conformal branes 013]. By deriving the appropriate boundary 
term for the action (|2. 1H . therefore, we automatically carry out the holographic renormalization of 
all these theories, whether asymptotically AdS or not, and without the need of any field redefinition. 
Explicit results for the case of non-conformal branes will be presented elsewhere |15j . 



2.1 Hamiltonian formulation of the variational problem 

The variational problem at infinity for the action (|2,ip is not well defined as it stands [8] [7j. For the 
case of asymptotically AdS gravity it was first shown in [16J that a certain asymptotic solution of 
the radial Hamilton-Jacobi equation renders the on-shell action finite. It was later shown [17] that 
the boundary term obtained by solving the radial Hamilton-Jacobi solution is the same as the one 
obtained via the standard method of holographic renormalization [18] [TH] [20], [21] [HI [22] [23] [24] ■ In 
[8] it was pointed out that the boundary term required to make the on-shell action finite is in fact 
the same boundary term required to render the variational problem at infinity well defined, while in 
[7] it was pointed out that this conclusion holds more generally, not just for non-asymptotically AdS 
gravity, but also for non-gravitational variational problems with a boundary at infinity, provided the 
variations are confined within a space of asymptotic solutions carrying a well defined symplectic 
form. Moreover, this boundary term always corresponds to a solution of the radial Hamilton- 
Jacobi equation. In practical terms this means that the leading solution of the Hamilton-Jacobi 
equation determined from the leading asymptotic form of the solutions must not contain transverse 
derivatives [7J. If this is not the case, it simply means that the space of asymptotic solutions 
corresponding to the chosen leading asymptotics is not well defined and it does not carry a suitable 
symplectic form. In such cases, one must first perform some kind of 'Kaluza-Klein reduction' to 
trivialize the transverse derivatives appearing in the leading radial asymptotics, derive an effective 
action for the KK fields, and then solve the radial Hamilton-Jacobi equation for this effective action 
of the KK fields [7] . The condition that the leading asymptotic form of the solutions should give rise 
to a solution of the Hamilton-Jacobi equation whose leading asymptotic form does not contain any 
transverse derivatives is automatically satisfied in the case of IHQCD, asymptotically AdS gravity, 
or non-conformal branes and we will assume that it is the case in the analysis below. 

To formulate the variational problem we start by writing the metric in an ADM decomposition 
[25], but with Hamiltonian time being replaced by the radial coordinate, r, emanating from the 
boundary at infinity. Namely, we write 

ds 2 = (N 2 + NiN l )dr 2 + 2N i drdx i + lijdafda? , (2.2) 

where N and JVj are respectively the lapse and shift functions, and 7^ is the induced metric on 
the hypersurfaces S r of constant radial coordinate r. The metric g^ u is therefore replaced in the 
Hamiltonian description by the three fields {N, ^,7^} on S r . In terms of these variables the Ricci 
scalar takes the form 

R[g] = R[j\ + K 2 - KijK ij + V M (-2Kn^ + n p V p n»), (2.3) 
4 See Appendix [B] for the explicit results in this case and [13] [14], where these results have been recently used. 
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where 12 [7] is the Ricci scalar of the induced metric 7^, the extrinsic curvature, Km, of the hyper- 
surface S r is given by 

Kij = ^ (7v " AiV, - DjN-) , (2.4) 

and -Dj is the covariant derivative w.r.t. the induced metric 7^. Moreover, K = j t3 Kij and 
= (l/iV, — iV l /iV), is the outward unit normal vector to S r . The total derivative term in this 
decomposition of the bulk Ricci scalar is an indication of the need for the Gibbons-Hawking term. 
Evaluating this term on T, r we see that it gives a contribution which is precisely canceled by the 
Gibbons-Hawking term. We therefore cirrivc ctt £i Lei grangian description of the dynamics of the 
induced fields {N, N{, 7^} on £ n namely 

2k 2 L = -j d d x^N (R[~i\ + K 2 - K)k( 

r ( 1 2N i 

+ ^ d d x^N{ — ^+Z(i P )x 2 )- 1 ^^d l ip + Z^)xd t x) 

+ (f* + ^r) + Z ^)diXdjX) ~ V{<p)\ . (2.5) 

Note that this Lagrangian involves no kinetic terms for the fields N and N{, which are therefore 
Lagrange multipliers, leading to constraints. The canonical momenta conjugate to 7^, (p and x are 
now obtained respectively as 

SL 1 

*x = 5 -^ = ^^m{x-N% X ), (2.6) 
while the momenta conjugate to N and iVj vanish identically. The Hamiltonian is given by 

H= f d d x (ir ij % + ir v <p + tt x x) - L = j d d x {NU + N-W) , (2.7) 

where 

W = 2.M(vrH-^ + I^ + V^(^) 

+^ (12[7] - dw&tp - ZMdix&X + V(<p)) , 
K = -2DjTT ij +TT^d i ip + Tr x d i x. (2.8) 

Hamilton's equations for the fields N and N{ lead respectively to the Hamiltonian and momentum 
constraints 

n = o, w = o. (2.9) 

Moreover, the symplectic form is given by 

Q = [ d d x (5TT ij A Sjij + 5tt v A 6<p + 5tt x A 5 X ) , (2.10) 



and it is independent of the value of the radial coordinate r [26J . 

The variational problem for the action (|2.ip can now be formulated in a regularized space, M To i 
whose boundary is defined as the surface S ro for some fixed r a . Provided r Q is sufficiently large, 
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T, ro is diffeomorphic to the boundary dA4 at infinity. Adding then a generic boundary term, Sb, to 
the action (|2.ip defined on A4 ro and considering a generic variation we obtain [7J 



S(S + S b ) 



d d+1 x(EOMs) + (L + S b 



Mr 



Sr 



+ 



d d x 



TT lJ + 



6Sb_ 

Hi 



Sjij + ( 7T<v, + 



5Sb 

Sip 



6<p+[ir x + 



5_Sb\ 

Sx) 



$X , (2-11) 



where the integrand of the integral over the bulk of M To is proportional to the equations of motion. 
In order for the variational problem at r Q — > oo to be well defined with boundary conditions imposed 
on the induced fields 7y, ip and %■> i- e - f° r a generic variation of the action with such boundary 
conditions to imply the equations of motion, one must demand that [7J 



(L + S b ) 



J +0. 



(2.12) 



The variational problem is then defined by considering variations of %j, ip and % within the space 
of generic asymptotic solutions of the equations of motion so that 



Sbl 



I To ' 



(2.13) 



where S is Hamilton's principal functional, i.e. a solution of the Hamilton-Jacobi equation where 
the values of the induced fields 7^ , ip and \ on S ro are totally arbitrary. It follows from (|2.1ip that 



7T 



'J 



5S 



hi 



TT, 



5S_ 

Sip 



7T 



X\r a 



ss 

~8x 



(2.14) 



where S is identified with the on-shell value of the action S on solutions with arbitrary boundary 
values for 7^, ip and % on ^r a - Including the boundary term S b in these relations leads to the 
canonically transformed momenta jTJ 



U lJ \ 



S(S + S b ) 



hi 



n 



s(s + s b ) 



dip 



n 



5{s + s b ) 



Sx 



(2.15) 



In order to determine the boundary term S b we need to determine the asymptotic form of 
Hamilton's principal functional. In other words we need to solve the Hamilton-Jacobi equation 
in a certain asymptotic sense, which we will specify below. The Hamilton-Jacobi equation can be 
derived from the following simple argument. Let S ro = S\ denote the on-shell action S with 
arbitrary boundary values for 7^, ip and x 011 ^r D - Then, 



dS To 
dr 

I27L41 dS ro 
dr 



+ 



+ 



iijh'Pix]-^ +^b>^x]^ +xh,v,x]j^ ) 

d d x {K lj % + n v (p + tt x x) 



+ H + L, 



(2.16) 



where d/dr Q denotes the partial derivative with respect to r Q . However, since S ro is the on-shell 
action, we must have S To = L and so we conclude that 



dS r 
dr 



+ H = 0. 



(2.17) 



For a generally covariant theory, like supergravity, the Hamiltonian vanishes identically since it is 
proportional to the constraints (|2.9p . It follows that, as a consequence of the general covariance 
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of the theory, the on-shell action does not depend explicitly on the radial coordinate, r, but only 
through the induced fields on the hypersurface S r . Moreover, we see that the Hamilton- Jacobi 
equation in a generally covariant theory is equivalent to the vanishing of the constraints, i.e. 



, i / / 1 \ 5S r 5S r 1 /SS r \ 2 1 i . . /<)W 2 ' 

= -^Vt(b[t] - SivS'v - z{'P)d iX a'x + v (f)) ■ 

The task of the next subsection will be to systematically solve these equations is a certain asymptotic 
sense. 

2.2 Recursive solution of the Hamilton-Jacobi equation 

We now turn to the task of solving the Hamilton-Jacobi equations (|2.18p in order to determine the 
boundary term Sb that makes the variational problem for the action (|2.ip at r a = oo well defined. 
As was argued in [7] and reiterated above, provided the variational problem is formulated within 
a well defined space of asymptotic solutions, the leading form of the boundary term, obtained 
as a solution of the Hamilton-Jacobi equation, will contain no transverse derivatives. It follows 
that the full solution of the Hamilton-Jacobi equation that corresponds to the boundary term we 
seek to determine admits an expansion in transverse derivatives. One can, therefore, try to solve 
the Hamilton-Jacobi equation by writing down an ansatz containing all possible terms allowed 
by general covariance at each order in derivatives. Although this approach, which was used in, 
for example, |T6] in the case of asymptotically AdS gravity, generically suffices for simple cases 
in low spacetime dimension and for a limited number of fields, it quickly becomes prohibitively 
inefficient and cumbersome. In particular, this approach not only unnecessarily includes terms in 
the ansatz that may happen to be absent in the particular theory, but also the number of equations 
one obtains for the undetermined functions in the ansatz is greater than the number of functions 
to be determined, and so many equations are redundant. Instead, the approach we will develop 
here is a systematic recursive procedure for solving the Hamilton-Jacobi equation, closer in spirit 
to the recursive method developed in [27] ■ In particular, the algorithm we will present computes 
systematically the nth term in the derivative expansion from lower terms, thus producing only the 
terms that do appear in the actual solution of the Hamilton-Jacobi equation. 

The basis of our algorithm is the following observation. Let us begin by writing Hamilton's 
principal function, S r , as 

S r = f d d x£( 7 ,v?,x). (2.19) 



From the relations (12,14j) then we have 

TT^djij + TT(pScp + tt x S X = 5C + dy{5-f, Sep, 6 X ), (2.20) 

for some vector field v l (5 , y, 5cp, 5x)- Since by construction the solution S r we are seeking admits a 
derivative expansion as r — > oo, S r increasingly approaches a solution of the form 



5(o) = \\ d d x^U(tp, X ), (2.21) 



for some function U((p, x)- Given the zero order solution (|2.2ip . we can compute corrections to this 
action in a systematic expansion in eigenfunctions of the operatoiH 

5 1 = jd d x2 7i ~, (2.22) 



3 In the case of asymptotically AdS or dS gravity coupled to matter, it was proposed in [T7] that one expands 
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namely (dropping the subscript r from now on) 



5 = 5( ) + 5(2) + 5(4) H , (2.23) 

where <5 7 5( 2n ) = (d— 2n)5( 2n ) • It is easy to see that the resulting expansion is a derivative expansion. 
Note that the operator (|2,22p agrees with the dilatation operator introduced in [T7] in the case of 
a constant potential, i.e. for the usual dilaton-axion system in conformal theories. However, the 
would-be dilatation operator for an arbitrary dilaton potential would lead to an operator whose 
eigenfunctions are highly non-trivial and so, from a practical point of view, would not serve as a 
good basis for expanding Hamilton's principal function. Moreover, an expansion in eigenfunctions 
of the operator ()2.22j) has the advantage of leading to algebraic in the induced metric jij equations 
for determining the terms 5(2 n ), which is the fundamental ingredient in our algorithm. 

Let us see how this works. Applying the general identity (|2.20p to the variation <5 7 we obtain 

2vr (2 „) = (d- 2n)C {2n) + d lV \ 2n) . (2.24) 

Since C is defined up to a total derivative, we can absorb the last term in this identity into £(2 n ) 
such that 

2vr (2 „) = (d-2n)£ (2n) . (2.25) 

The significance of this relation will become clear shortly, when we write down the equation de- 
termining £(2n)- Since f|2.25|) holds only for a certain choice of total derivative terms in £(2 n )> in 
principle we should keep track of total derivative terms as well. However, as we shall see, this 
will not be necessary in our recursive procedure. In particular, we will determine C^n) a * each 
order up to total derivative terms. The canonical momenta at this order can then be obtained 
by differentiating 5( 2n ) = / d d xC( 2n y Since total derivative terms do not influence the canonical 
momenta, we will still get the correct expressions for the canonical momenta. 

Inserting the leading term ()2.21[) of the above expansion in to the Hamilton-Jacobi equation 
(|2.18p we find that the function U(<p, x) satisfies the equation 

{d^Uf + Z-\y){d x Uf - -^-U 2 + V(<p) = 0. (2.26) 

However, using the leading term (|2.2ip in the relations (|2.14p and identifying the momenta from 
the expressions (|2.6|) we find that to leading order the induced metric takes the form[f| 

Jij = e 2A g m {x), (2.27) 

where g(o)iji x ) * s an arbitrary metric on the boundary and the fields A, ip and x satisfy the flow 
equations (|1.7p . but with W replaced with U. But then, these first order equations in combination 
with (|2.26p imply the second order equations (jl.6p . As we have seen above, the second order 
equation for the axion is integrable, giving in this case 

= &-«", (2.28) 

dx 



S r and the corresponding canonical momenta in eigenfunctions of the dilatation operator, Sd, obtained from the 
relation d r = J d d x (jijj^— + Ylffjjj r ^°°> jSd, where / denotes generic matter fields. In that case expanding in 
eigenfunctions of the dilatation operator is the most efficient recursive method to solve the Hamilton-Jacobi equation, 
because it amounts to solving the zero order problem (|2.26l) and the one determining the higher derivative terms 
simultaneously. However, in more general cases the operator corresponding to the leading asymptotic behavior of 
d r is not very useful in practice since its eigenfunctions are not simple functions of curvature invariants and matter 
fields. In most cases it is easier in practice to first solve the zero order problem (|2.26[) to determine U(<p, x) an d then 
to expand in eigenfunctions of <5 7 . Of course in the case of asymptotically AdS or dS gravity the two approaches 
produce identical results, even though the two expansions differ order by order, the difference being that the latter 
expansion resums all zero derivative terms into the leading term. 
6 We gauge-fix the lapse and shift functions to N = 1 and Ni = 0. 
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which is the analogue of (jl.lOp , We therefore see that any \ dependence in U leads to a finite 
contribution in Hamilton's principal function S and hence we can set the integration constant 
c = 0. We conclude that the function U (p) required to make the variational problem well defined 
can be taken, without loss of generality, to be independent of the axion \- As we will see shortly, 
this leads to a significant simplification of the analysis to determine the required higher derivative 
terms in S. 

In order to determine the leading solution (|2.2ip of the Hamilton-Jacobi equation, therefore, 
one needs to solve the reduced equation 

(d v U) 2 - 1 ^ T U 2 + V(ip) = 0. (2.29) 

For an arbitrary potential this equation can be transformed into an Abel's equation of the first kind 
[IT] . which is generically non-integrable. However, we need not find the general solution of this 
equation. Any solution of this equation that ensures that cp has the desired general asymptotics 
via the relation 

<p = ±^ 1 ± T U^^-~V^p), (2.30) 

suffices. In particular, if U is such a solution and it is not isolated in the space of solution^ 
then it is easy to see that any solution in the vicinity of U is of the form U + eAU, where e is an 
infinitesimal parameter and AU = 0(exp(— dA)). Hence, the difference between two such solutions 
only contributes a finite term in 5(o)> an d it is therefore irrelevant. The reason for restricting this 
argument to the vicinity of the original solution U is that there exist solutions at infinite parametric 
distance from U that change the leading asymptotic behavior of (p. Such solutions are excluded by 
the requirement that ip has the correct asymptotics. 

Now that we have determined that U((p) is independent of the axion, inserting the above 
expansion of Hamilton's principal function in the Hamiltonian constraint and matching terms of 
equal 5 7 eigenvalue we obtain for the higher order terms 

U '^li I ^ xC W> ~ (^Zt) C/ ^) £ (2-) = n <?*)> n > 0, (2.31) 



where 



1 



K {2) = -^^(R^l-di^cp-Z^dixd'x), 

n—l / 

~2 \ I ■77V„.__\*.'n-/v./.. ...w? — 

d-1 



i " _1 / -1 

K(2n) = -2K 2 7 _ 2 I VT( 2m )j-Vr(2( n _ m ))i - - _ 1 7T(2 m )7r(2( Tt _ m )) 



m=l 



"7 7r V(2m) 7r ^(2(n-m)) + ~{ Z 1 (f^X (2m) 7r x(2(n-m)) ) > U > 1. (2.32) 



Importantly, these are linear equations and only involve a derivative w.r.t. the dilaton, <p, and not 
the induced metric 7^ or the axion. The absence of a derivative w.r.t. the induced metric is due 
to the relation (|2.25p . while the absence of a derivative w.r.t. the axion is because we have shown 
that U only depends on the dilaton. We will now solve these equations up to the order required in 
order to determine the boundary term Sb for d = 4. In the analysis below we will take U'{(p) ^ 0, 
which is guaranteed by the requirement of a logarithmic running for the dilaton. However, the case 
of a constant potential V and, consequently, constant U arises in conformal theories such as J\f = 4 
super Yang-Mills and so it is interesting on its own right. For completeness we discuss this case in 
Appendix [Bj 



7 See Appendix [B] for an example where this happens. 
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Turning now to the case U'(tp) ^ 0, we notice that the linear equation (|2.3ip for Cr 2n \, n > 0, 
admits the homogeneous solution 

4^ = ^n)[7 ! x]e XP ((^t) fj^ff^) ■ ( 2 - 33 ) 

where ^ r (2n)[7)X] is a covariant function of the induced metric and the axion of weight d — In. 
However, this solution contributes only to finite local counterterms and can be ignored. To see this 
notice that the leading order solution (I2.2ip implies via the Hamilton-Jacobi relations (j2.14j) that 
the induced metric takes to leading order the form (|2.27p with 

Now, by construction, J-(2 n ) [7> x] nas weight d — 2n and so 

C (2n) ~ e {d - 2n ^ A x e -("M = finite. (2.35) 
We are, therefore, only interested in the inhomogeneous solution of (|2.31|) . In other words, formally, 

C (2n) = e-(")^) f J^ e ^ A ^n {2n) m. (2.36) 

This integral is well defined if TZ^n) does not involve derivatives of the dilaton, ip, but it requires 
some caution when it does. In general, we can write more precisely 



£(2n)=e^ d - 2n ^F {2n) , 



(2.37) 



where F^n) satisfies 



(2.38) 



for some vector field v r 2n \ l ((p, Sip). In Table [U we have listed the local functional F^ 2n ^((p) for a 
number of generic source terms 7Z( 2n j (ip) that we will need for our calculation. A detailed derivation 
of the formulas given in Table [1] is provided for the reader's convenience in Appendix El Both in 
Table [1] and in Appendix [A] we make extensive use of the short-hand notation 

L 9 _^m g _ (d _ 2n)A I"? 'd£ e (d-2n)A ^ym ^ ^ 
Jn,m J U 

where A(<p) is given by (|2.34p . 

The integration formula (|2.38|) allows us to develop an algorithmic procedure for evaluating 
Hamilton's principal function, S, iteratively. Namely, given the source term of equation (|2.3ip at 
order n, (I2.38P is used to obtain Ct 2n \ up to an irrelevant total derivative term. Differentiating this 
with respect to the various induced fields gives the corresponding momenta at that order. Finally, 
using these momenta, as well as those of lower orders, one determines the source term at order n + 1 
via (|2.32p . The procedure is then repeated to the desired order. This algorithm is schematically 
outlined in Fig. [TJ Let us now carry out this general algorithm to order n = 2, which is sufficient 
for evaluating the boundary term for the action (|2.ip in d = 4. In general, at each order n the 
source term, TZ( 2n ) and the corresponding inhomogeneous solution, £( 2n ), of the functional equation 
(|2,3ip can be written respectively in the form 

1 N n 

K (2n) = -^v^E^) 7 ^' ( 2 - 4 °) 
i=i 
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e - (d - 2n) A F(2n) 


r lm (<p)t ili2 - im d il (pdi 2 <p...d im (p 




r 2 (ip)t i W i D j i P 




(r l22 ((p)t^ kl + s l22 (<p)4 J ' W J dnpdjipD k Diip 


fn,3 s 1 2 2 (ip)t t J kl d i ipd j ipD k Diip 


(r 2 i(tp)t^ kl + s 2 2{(p)t 2 jkl ^ DiDjtpD k Dnp 


-2f£ 3 U>A% (£) f n ^)q k d^D k D^ 



Table 1: The result of the functional integration described by (|2.38p for various source terms. 
Here, t %lt2 '" %m and t %1 are arbitrary totally symmetric tensors independent of ip, while = 
I( 7 *V 7 + 7 tf y fc + 7 ii 7*')> l T l = \ [l ik l jl + 7^7 J ' fc - 2 7 ^7 fc/ ). Moreover, F {2n) is given up to 
terms of the form e^ d ~ 2n ^ A DiU l , for some vector field W . Since Ci 2n \ is related to F( 2n ) as i n (|2.3T|) . 
W corresponds to a total derivative term in £(2n) and so we need not determine this term explic- 
itly. Note that the source term r l 2 2 only contributes to W and so it can be ignored. A detailed 
derivation of these results is given in Appendix lAl 



n 



(2n) 



C 



(2n) 



{ 7r (2n)} 



{7T 



(2n+2) 



} <" 



(2n+2) 



ft 



(2n+2) 



n 



(2n+4) 



(2n+4) 



Figure 1: A schematic illustration of the algorithm used to determine Hamilton's principal function 
iteratively. Here, the operators j and 5 stand respectively for the functional integration defined by 
the formula in (|2.38|) and for functional differentiation with respect to the induced fields. 
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I 


Tl 


4(<p) 




1 


R 


i 




2 




-i 


f£(-l) = -MM 


3 


dixd l x 


-Z{fp) 


JJ (-z(v)) = -eM 



Table 2: Summary of the source terms (|2.40p and the corresponding solution (I2.41D of the Hamilton- 
Jacobi equation at order n = 1. The source term is determined iteratively as in (|2.32j) . while C^n) 
are determined by solving the linear equations (|2.31|) via the integration procedure described in the 
text. 



and 

, N n 

ZK 1=1 

where 0^(92) and V^i^p) are scalar functions of ip and are quantities involving fields other than 
(p, as well as derivatives of (p. The number N n is the number of such quantities at each order n. 
For n = 1 the source 7£(2) is given in (|2.32|) . Applying the results in Table [1] to this sources we 
obtain £(2)- The result is summarized in Table El 

In order to move to the next order we first need to evaluate the canonical momenta at order 
n = 1 by differentiating £( 2 ) with respect ot the induced fields. After a little algebra we obtain: 

7r (2) tf = -i,/7(W -E'D'D^+^M -2S")d i ^i P + ^ed i x& i x 

^( 2 ) = -^v^QM'^V + Mn^-ie'^x-s'i? 

vr x(2) = -^^(en 7X + G'^x). (2.42) 

Inserting these in the expression for 7^(4) in (|2.32p results in an explicit expression for the source 
term at order n = 2. This is a rather complicated source, involving 20 different terms, and can be 
read out from the first three columns of Table [3j Using the results of Table [1] we then obtain the 
fourth column in Table [3l giving £(4). This table is the main result of this paper providing together 
with Table [21 the boundary term necessary to make the variational problem of a generic action of 
the form (|2.ip well defined in any dimension up to and including d + 1 = 5. 

Before we apply this general result to IHQCD, a technical comment is due. The boundary terms 
listed in Tables [2] and [3] generically have poles at d = 2 or d = 4. This happens whenever there 
is a conformal anomaly, or, in terms of the bulk language, whenever the boundary term required 
to make the variational problem well defined breaks the bulk diffeomorphisms corresponding to 
translations in the radial coordinate. The way to handle these poles is to relate the radial cut-off, 
ro, to the parameter d, a^lro = l/(d — d*), where d* is the dimension of the boundary. An example 
of this issue arises in the case of a constant dilaton potential, discussed in Appendix IBl 

8 The radial cut-off ro and the boundary dimension d are arbitrary parameters in the radial Hamiltonian formalism. 
However, the fact that this formalism, without relying on any additional input from the asymptotic expansions, can 
handle conformal anomalies in d* dimensions by relating the radial cut-off to the parameter d suggests that a radial 
cut-off regularization corresponds to dimensional regularization in the dual field theory. Taken at face value, this 
identification has profound consequences for the physical significance of the holographic direction. 
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T 1 

>2 




^(*0 


1 




45 2 


■f v 45 2 


2 


R 2 


^/2 d ^2 


fV /"^/2 d ^2\ 
J2,(A" d-1" / 


3 


D i D^<pD i Djip 


4S' 2 


4/2 V 2 -' 2 


4 




-4S' (M - 2H") 


| g 3 (45'(35" - M) - MM' - 2^7',4'd 2 , (^) f| 2 (65' 2 - M 2 )) 


5 




MM' + 25' (M - 45") 


-1 f 2 3 ( 4S ( 3S - M ) - MM - 2C/ M 5 I (w) f 2 2 ( 6S - M )) 


6 




M 2 - 45' 2 


£(M 2 -45' 2 ) 


7 




i M /2 +^jM 2 2M5" 


fi(|M' 2 + " 2M5") 


8 




| (3^-4)^2 




9 






-8 £55' 


10 


R ij diipdjip 


45 (M - 25") 


4/& (5(M - 25") + 2U'A% (£) f£ 55') 



11 




-23' (M - 23) 


-2g 1 3'(M-23) 




12 


Rdiipd z ip 


-S'M' - j^SM + 433" 


- E'M' - ^3M + 433" + 2U'A% (£) 3'(M - 


23)) 


13 


R ij dixdjx 


43G 






14 


D i D j tpd i xdjX 


-43'0 


-4/^3'Q 




15 


(d^x) 2 


z- 1 e /2 + 2(M- 23") e 


^ 2 (Z- 1 /a + 2 (M — 23") 9 + WA'dl (£) £ S'Q) 




16 


U 1 ipd i xd t x 


23'e - Me' 


(23'G - M6') 




1 7 
1 1 


Rdixd l x 


- U - 


J2,o^ ° ~ a=T" u ^ 




18 


diipd l (pdjxd j x 


x e (43" - ^M) - |M'G' 


I f£(G(43" - gzjM) - M'G' - 2A'l7'fl£ (£) f£(23'G - 


Me')) 


19 


(□ 7 x) 2 








20 


Ojxdifd'x 


2z- 1 eo' 


£ 2Z-1GQ' 





Table 3: This table summarizes the source terms (I2,40p and the corresponding solution (|2.41j) of the Hamilton- Jacobi equation at order n = 2. The 
solution Cu\ of the Hamilton- Jacobi equation at order n = 2 can simply be read off the fourth column via the relation (|2.4ip . This table is the 
main result of this paper providing, together with Table [21 the boundary term necessary to make the variational problem of a generic action of the 
form (|2.ip well defined in any dimension up to and including d + 1 = 5. 



3 Application to IHQCD 



The results presented in the previous section are applicable to any action of the form (|2,ip . In 
particular, Tables [5] and provide the general boundary term that renders the variational problem of 
any action of the form (|2.ip well posed for boundary dimension up to and including d = 4. However, 
additional assumptions that enter into specific models often lead to significant simplification of this 
boundary term. In this section we will apply these general results to the specific problem of IHQCD 
|10j described in the Introduction. The significance of this boundary term will be demonstrated by 
deriving asymptotic expansions analogous to the Fefferman-Graham expansion for asymptotically 
AdS gravity [6], providing general expressions for the one-point functions of the dual operators 
in terms of the coefficients of these expansions, and finally giving the correct holographic Ward 
identities for IHQCD. 

In deriving the boundary term in Tables [2] and [3] we have assumed that the leading asymptotics 
of the induced fields follows from a leading solution (|2,2ip of the Hamilton-Jacobi equation which 
does not involve transverse derivatives. However, no specific asymptotics was assumed. Taking 
into account the leading asymptotic form of the induced fields of the IHQCD model allows for 
various simplifications of the result of the previous section. To see what simplifications occur let 
us start by writing ip = £logA, for some non-zero constant £. Then, by abuse of notation, we will 
correspondingly write U(X), A(X), etc. as functions of A. As we have argued above, U is a function 
of the dilaton A only and not of the axion. Moreover, it is a solution of equation ()2.26|) . which now 
becomes 

_4_^_ r ^Q 2 . m . (3 ,) 

At this point we need to invoke some information about the asymptotics. First, recall that we 
want to consider potentials of the form (|1.4|) as A — > 0. Moreover, U(X) defines the leading form of 
the metric and dilaton asymptotics via the flow equations (these are the same as (jl.7p but with W 
replaced by U, and follow from combining the relations (12. 6p and (I2.14p for the canonical momenta) 

The leading metric and dilaton asymptotics adopted in the model of IHQCD is |10j 

A^T 1 , A~ -bol^X 2 + b 1 r 1 X 3 , as A^O (3.3) 

where the constants &o and b± were argued in |10] to be related respectively to the one and two- 
loop perturbative beta function coefficients. According to that argument, subleading corrections to 
the asymptotic form of the dilaton correspond to higher loop coefficients in the perturbative beta 
function, which are scheme dependent. In that sense, the terms in (13. 3[) are universal in the model 
of IHQCD, and any dilaton potential V(X) chosen must be such that it is compatible with this 
leading asymptotics in the UV, while subleading terms in the UV expansion of the potential can 
differ for different choices of the dilaton potential. 

One should keep in mind, however, that the beta function of any operator whose coupling 
transforms inhomogeneously under scale transformations is scheme dependent and, therefore, not 
a physical observable. In other words, the beta function of any operator that acquires anomalous 
dimension is scheme dependent. For such operators the scaling dimension 70 is not constant along 
the RG flow. Only the value of this scaling dimension at the fixed points is a physical observable. 
Nevertheless, the combination /?£> • 0(x) is a physical observabl^l, since this combination appears 
in the trace of the stress tensor, whose scaling dimension does not renormalize. As we shall see 
below from the trace Ward identity, the bulk dilaton A should be thought of as the holographic 
dual to the operator /3(Xym)^ F 2 . Isolating the beta function factor in this operator is a scheme- 
dependent procedure, but the product transforms homogeneously under the renormalization group 

9 We are grateful to Elias Kiritsis for pointing this out to us. 
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flow. It is precisely this combination that appears in the trace Ward identity. In fact, this is a 
generic mechanism of how one can accommodate operators with running scaling dimensions in a 
supergravity setting, thus going beyond the realm of operators with protected dimensions. 

In order to ensure that the leading asymptotic behavior of the metric and dilaton is of the form 
(|3.3|) we take the function U(X) to be of the form 

U(X) = -^-^X + - i U 2 X 2 + - i U a X a + o(X a ), a>2, (3.4) 

where 61 = 2£~ 2 U2, and we have introduced a term of order X a , a > 2, in order to allow for 
potentials with non-integer powers of A, as were considered e.g. in 128] . However, it should 
be emphasized that the full closed form of the function U (A) is required in order to construct 
the necessary boundary term that makes the variational problem well defined. In particular, to 
write down explicitly the boundary term one necessarily needs an exact solution of (|3. 1 H . Such 
an exact solution is of course dependent on the particular choice of dilaton potential. Since we 
want to keep the discussion general here, we will not fix the dilaton potential or the function U(X), 
beyond the requirement that as A — > it behaves as in (|3.4p . Once a dilaton potential is specified, 
one then will simply need to find an exact solution of (I3.ip that behaves asymptotically as in (I3.4p 
and use it in the expressions below, all of which are expressed in terms of a generic U(X). Of 
course, such a solution will not be unique since it can be shown [29] that if Uo(X) is such a solution, 
then there is a continuous family of deformations of this solution, with the deformation behaving 
as X du exp(— d/boX) as A — > 0, where the value of the exponent v is given below. However, any 
member of this continuous family of solutions is equally good for the purposes of evaluating the 
boundary term since any difference arising from this deformation only contributes finite terms to 
the boundary term. 



3.1 Boundary term for IHQCD 

In this section we will make use of the asymptotic form (|3,4h of the function U(X) in order to 
simplify the general boundary term derived in the previous section. To this end we start by noting 
that, depending on the value of the exponent a, the function ^4(A) defined in (I2.34D has the following 
asymptotic behavior as A — > 



X~ u e b o x 1 



A-W[l-i 



\->e*» I 1 - £ 



aUa 
2U 2 N 2 



X a ~ 2 + (A 



a-2 



^ + gMA + o(A) 



2u 2 y 



a < 3, 



a 



u 2 



UA + o(A) , a>3, 



(3.5) 



where v = -r-r + Wii- Using these asymptotic expansions and the integral identity 

" 1 k "0 



f 



dX'X'^e^ 



h^^ fe^o ^ (^) n A" + O (A-+ 1 )) , u > 0, Vn G e N, 



=0 r<» 
-^-A^" 1 + const., 
log A + const., 



u = 0, u^l, 
uj = 0, u = 1, 



(3.6) 



we can determine the asymptotic form of all functions listed in Tables [2] and [3l These functions 
involve integrals of the form 



X^^f (XA'(X)) m c^ d ~ 2n ^ 



dX 



X 2 U'(X) 



=("") A W(Ai'(A))- m A A (3.7) 



whose asymptotic form we tabulate in Table [H With these results one can now determine the 
asymptotic form of the functions in Tables [2] and [3l which we present respectively in Tables [5] 
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d-2n 


A + m 


C x \A 

Jn,m 


> 


any 


^(i+(I^-£t)m+^)) 





7^0,1 


-Mf^y (l + ^ - 1 + 2m)|| + A + 0(A')) 





1 










^(l-(f^+A6 ,)AlogA + O(A)) 



Table 4: The asymptotic form of the integral (|3.7p for a function £7 (A) that behaves as in (|3.4j) as 
A -> 0. 

1 R -2H(A) = / I (l-g^A + 0(A 2 )) 

2 e 2 A- 2 a t A^A -M(A) = -^(l+(^-^)6 A + O(A 2 )) 

3 da** -Q(X) = - 1 ^(M^y 1 (z +[z 1 -^Z )x + O(X 2 )) 

Table 5: The leading form of the solution (|2.4ip of the Hamilton-Jacobi equation at order n = 1 for 
IHQCD. The third column gives the asymptotic form of the terms in the fourth column of Table [2j 

and [6l In deriving these asymptotic forms we have made repeated use of the crucial fact that 
diX = 0(A 2 ). This follows from the general asymptotic form of the dilaton A and will be derived 
in the next section. 

Note that even though we have listed the leading asymptotic behavior of the order n = 1 
terms in Tabled we actually need the exact closed form expressions of the functions M(A), H(A) 
and 0(A) in the boundary term. The reason is quite obvious. Namely, the expansion of the 
solution of the Hamilton-Jacobi equation in eigenfunctions of the operator <5 7 , i.e. the derivative 
expansion, is an expansion in powers of e~ 2A ~ X 2u exp(— 2/boX). This expansion, therefore, is 
non-perturbative in A. However, at each order in the derivative expansion the boundary term 
contains an infinite expansion in powers of A. In other words, the expansion of the boundary term 
is a double expansion! 10 ! Clearly, at orders n = and n = 1 in the derivative expansion we must 
keep the entire perturbative expansions in A, since there are divergences coming from any power of 
A no matter how large. This is the reason why the full closed form expressions for U{X) at order 
n = and of M(A), S(A) and 0(A) at order n = 1 must be kept. At order n = 2 in the derivative 
expansion, however, this is no longer necessary. Indeed, for d = 4 the factor e 4j4 coming from the 
volume element ^fy exactly cancels the factor e~ 4A coming from the four-derivative terms. Hence, 
the divergences at order n = 2 are only power- like or logarithmic in A. The first few terms of the 
asymptotic form of the boundary terms listed in Table EJ therefore, suffice. 



1 Generically, there are logarithmic in A terms as well, but these can be included with the powers in A since they 
do not affect the counting of the derivative expansion. 



19 



T 1 

'2 
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3 {X^DW^XDiDjX - 2\^d i \d j \D i D'J\ + A- 4 («9;Ad*A) 2 ) 

4 e 3 {X^diXdjXDW^X - A" 4 (9iA^A) 2 ) 
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ZoX- 1 + 



2g 2 bp _ 2U 2 
d-1 ^ 



Zo-ZOlogA + ^A ) 



^(JMg,JVg)-l 
(d-2) 2 (d-l) Z + <-AA) 
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(d-2) 2 ft 
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4^Z (A- 1 + O(A°)) 



^2-§t)^0 + O(A) 



H^oA-+ (m-^ )Z - Zl logA + O 



2(d-2) 2 (d-l)fe 



d£3 (M 3^2)-l 
"2(d-2) 2 (d-l)feo Z ° ^ A + U( ^JJ 
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2g 2 bo _ 2C/ 2 
d-1 Fto 



Z - Zi ) log A + O 
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Table 6: The leading form of the solution (|2.4ip of the Hamilton-Jacobi equation at order n = 2 and d = 4 for IHQCD. The third column gives the 
asymptotic form of the terms in the fourth column of Table O 



The entire boundary term that renders the variational problem for IHQCD well defined, there- 
fore, takes the form 



glHQCD _ _ + s ^ + s ^ 



+\a(X) (R t ^ - l -R 2 ^ + \b(X) - d iX d jX - \(U lX ? 



(3.8) 



where we have set d = 4 and introduced functions a(A), 6(A) and c(A), which are given by 
a(A) = ~ (6 - 1 A- 1 + 2ClogA + c 1 ), 

6(A) = PiMffiyiZo (b^X- 1 + (2C - log A + c 2 ) , 

c(A) = -|(M p 3 iV c 2 )- 2 Z 2 ^A" 1 + (2C - log A + c 3 ) . (3.9) 

Here 

£ 2 C/ 2 

and c\, C2 and C3 are arbitrary constants. The terms they multiply correspond to finite local 
counterterms and they reflect the usual scheme dependence. The functional derivative of each of 
these three terms w.r.t the induced metric gives a local transverse and traceless tensor. These 
three tensors are given in ()3.33[) - (j3.35[) below. Moreover, the tilde in the n = 2 solution, Su\, of the 
Hamilton-Jacobi equation is there to remind us that this is not the full solution of the Hamilton- 
Jacobi equation at order n = 2. Namely, «S( 4 ) contains only the local divergent part (plus scheme 
dependence) of the n = 2 solution. In addition, there is a finite part, 5(4), that corresponds to the 
renormalized action and which we have not determined. Recall that all the homogeneous solutions 
of the equations in the previous section were ignored precisely because they contribute only to 
this finite part, which cannot be determined from the derivative expansion of the Hamilton-Jacobi 
equation alone. As usual, to determine this part one must impose some regularity or boundary 
condition in the deep interior of the spacetime. 

The form of this boundary term deserves some close inspection. By far the most remarkable 
feature is that it is totally covariant. In particular, contrary to the boundary term (|B.12|) for the 
strictly asymptotically locally AdS dilaton-axion system, it does not break the bulk diffeomorphisms 
corresponding to shifts of the radial coordinate. A direct consequence of this fact is that, as we 
shall show below, the trace of the stress tensor does not depend explicitly on the sources, but only 
implicitly through the dilaton one-point function. In other words, there is no conformal anomaly 
in IHQCD, although conformal invariance is broken via the dilaton one-point function. A related 
observation is that there are terms proportional to A -1 ~ r at order n = 2 in the derivative 
expansion. These terms are in fact nothing but the conformal anomaly for a gravity-axion system 
in asymptotically AdS space, as can be seen from (|B,12|) in Appendix [Bj This picture suggests 
that the field A effectively acts as a compensator for scale transformations or shifts of the radial 
coordinat qjj. Later on we will confirm this by showing that the source of the dilaton A can be 
removed by a bulk diffeomorphism that induces a Weyl rescaling of the boundary metric, while the 



n We are grateful to Kostas Skenderis for pointing out this interpretation to us. 
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one-point function of the dilaton contains exactly the would-be conformal anomaly of the gravity- 
axion system. 



3.2 One-point functions and Fefferman-Graham expansions 

Having solved the Hamilton-Jacobi equation for IHQCD asymptotically, we are in a position to 
derive the full asymptotic behavior of the bulk fields, i.e. the generalized Fefferman-Graham ex- 
pansions, without solving asymptotically the second order equations of motion. Not only is this 
approach considerably more efficient than directly solving the second order equations of motion, 
but it also avoids the necessity of making an ansatz for the form of the asymptotic expansions. As 
we shall see, the structure of these asymptotic expansions can be highly involved and practically 
impossible to guess a priori. However, the asymptotic solution of the Hamilton-Jacobi equation that 
we have obtained above already contains all information about the form of the asymptotic expan- 
sions. Additionally, it will automatically tell us how the one-point functions, i.e. the renormalized 
momenta, are related to the coefficients of the asymptotic expansions. 

The key ingredient in deriving the asymptotic expansions from the asymptotic solution of the 
Hamilton-Jacobi equation is the first order flow equations 

i = K i Z -> W ±^, (3.11) 

which follow from combining the two expressions for the canonical momenta given in (12. 6p and 
(|2.14p . Now, the solution of the Hamilton-Jacobi equation we have determined above takes the 
form 

5 = 5(0) + 5(2) + 5(4) + 5(4), (3.12) 

where 5(4) is undetermined and remains finite when the radial cut-off is removed, while the other 
terms are given in f)3.8f) . Inserting this form of the solution in (|3.1ip we obtain the following flow 
equations. 

A = C 2 >?d\U 

+-A" 1 (2M - Xd x M) di\d i \ - MD 7 A + -^^dx^d^x + C 2 X 2 d x ^R 



■2 " ' ' 2 

+ dxc&xd'x) 2 



+7^ A 2 (dxa{R i3 R lj - \r 2 ) + d x b 



4£ 



( Rlj - l -R ll3 )d\dix-\(^X? 



+K 2 r 2 A 2 -^7fA( 4) , (3-13) 



X 



-z- 1 (on^x + dxQdiXd'x) 



I 2/7—1 1 ^ 



^Rrafi&x + ~A(ra 7X ) + 2cd kX d k xd lX 



(3.14) 
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2 

-AER^ + AdxEDiDjX + 4 (d 2 x E - ^ 2 A~ 2 m) diXdjX - 2Qd iX d jX 

(sR + ^M\~ 2 ed k \d k \ + l -Qd kX d k X 
-2aR ikjl R kl + ^aRRij + 2D {l D k {aR k j} ) - 0,(0%) - ^D^aR) 
{^a{R kl R kl - ^R 2 ) - ^D k Di(aR kl ) + ^D>i?)) 

+ \bli 3 {Rki ~ \Rlki)d k X d l X ~ bR {lk d k X d 3)X + \bRdixd.X ~ bR ikj id k X d l x 
+ l -bR i] d kX d k x + D {l D k (bd kX d j)X ) - ^(bdtxdjx) ~ \-yijD k D l (bd k x&x) 
~D i D j (bd kX d k X ) + ^i j 1 (bd kX d k X ) + bDiD jX n iX ~ D {l {bU lX D j)X ) 
+^ lij D k {bU lX D k x) ~ ^ijK^x) 2 + c {^H j d k xd k X ~ 2d lX d jX ^j d lX d l X 
+4k 2 yjikjji - ^pklli^j ^=%) M . (3.15) 

Here, 7r\(4)i ^x(4) an( ^ ^(i) M denote respectively the functional derivatives of Su^ w.r.t. the dilaton, 
the axion and the induced metric. Since 5(4) is undetermined, so are these renormalized momenta, 
which define the one-point functions of the dual operators [17] . The presence of these terms in the 
flow equation will lead directly to the identification of the normalizable modes in the asymptotic 
expansions. 

Although these flow equations look rather complicated, they can actually be solved in a fairly 
straightforward way by noticing that the asymptotic solutions will in fact be two-scale expansions. 
Namely, there is an expansion in exponentials of r coming from the derivative expansion. Moreover, 
at each order in this expansion there is an expansion in powers and possibly logarithms of r. Since 
the structure of the expansion in exponentials of r is pretty clear from the derivative expansion, it 
is useful to separate the two expansions by writing explicitly 

1ij(r,x) = e 2r/i (p/( )ij( r > x ) + e ~ 2r/i 7(2)ij( r , x ) + e _4r ^7(4)y 0*, x) H ), 

A(r, x) = A (0 ) (r, x) + e~ 2r/i X {2) (r, x) + e~ 4r/e A (4) (r, x) + ■ ■ • , 

X (r,x) = X(o)i r ,x) + e' 2r/e x(2){r,x) + e~ 4r/e X(4) {r, x) H , (3.16) 

where the coefficients of the exponentials here are undetermined functions of r, only constrained 
by the requirement that their asymptotic expansions only contain powers and logarithms of r, but 
not exponentials. 

Inserting these expansions into the flow equations leads to non-linear first order equations for 
the order zero coefficients and linear first order equations for the higher order coefficients. Namely, 
for the induced metric we get 

7(0)ij + 2 (j + ^(A( ))J 7(0)« = 0, (3.17) 
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2 , 

7(2)*j + (A( ))7(2)*j = - 4S (A(o))^[7(o)] + (A (0 ))-D(o)i-D(o)jA (0 ) - 2e(A (0 ))<9iX(o)djX(o) 



+4 (s"(A (0) ) - ^ 2 A (0 )- 2 M(A (0) )) ^A (0) ^A (0) 

+§T(o)v (-f/ / (A (0 ))A (2) + S(A (0) )i2[7(o)] (3-18) 

1 £ 2 

+ 9©(A(o))^ ) (o)fcX(o)^ > (o) X(o) + irA(o)~ M(A( ))-D( ) A( )£>(o)feA( ) 



/2 2\ 2 / 1 \ 2 

7(4)ij + ( ^ - j J 7(4)« = - g ( A(4)^' + 2 A (2) 2f/// ) 7(0)ij - gA (2 )[/ / 7(2)ij 

-4A( 2 )H / i?i i [7( 0) ] - 4~ f J D( )fe-D(o)(i7(2)f) - 2 D (o)7(2)y - 2^(0)^(0)^7(2)^ 

+4A(2)H / 'L>( ) i L>( ) i A(o) + 45' ^D( )iD(o)jA( 2 ) + Q^o)*^)*? ~ - D (o)(i7(2)fei)^) £> (o) fc A (0 )^ 

+2 - yA (0) - 2 i\/) %A (0 )^)A (2) + 4A (2) (V - y\o)" 2 (M' - 2A (0) - 1 M)) ^A (0) ^A (0) 
-2\ 2 )®' diX{o)djX(p) - ©%X(o)<9j)X(2) 

+ g7(2)ij ( Si2[7( )] + yA (0) ~ 2 ML> (0)fe A ( o)£>(o) fe A (0 ) + ^QD {0)k X(o)D (0) k X(o) 
2 



2 / \ c 

+ g7(o)*i A (2 )H'i2[7 ( o)] + H [D [0) k D {0) l 1{2)kl - 0(0)7(2) - 7(2) fc '-R[7(o)]fci) + yA^^x 
(a (2 )(M' - 2A (0) - 1 M)Z) (0)fe A (0) Z) (0) fc A (0) + 2MZ) (0)fc A (0) l} (0) fc A (2) - M 7(2) fci a fe A (0) ^A ( o)) 



+\\2)®'D { Q )k X{o)D(o) k X{o) + ®D (0)k X(o)D {0 ) k X(2) ~ (2) kl d k X(o)diX{o) 

2 

-2aR ikj i [7 ( o) ]i? fei [7(0)1 + gO-R[7(o)]-Rij[7(o)] + 2#(o)(^(o)fc(a^ fe j)[7(o)D - □(o)(ai^[7 (0 )]) 
-|^(o)<^(o)j(a^[7(o)]) + 7(o)« Qa(^fci[7(o)]^ fc '[7(o)] " £# 2 [7(o)]) - \D {0)k D {0)l (aR kl [j {0) }) 
+^ n (0)(a-R[7(0)])^) + ^&7(o)y ^-Rfci[7(0)] - ^[7(o)]7(o)fci^ ^(o) fc X(o)^(o)'x(o) 
-bR( ik h(o)} D (o) k X(o) D (o)j)X(o) + ^ b R[j(o)}diX(o)djX(o) ~ bRi kj ib( )} D (o) k X(o) D (o)X(o) 
+^ b Rij[l(o)} D (o)kX(o)D( ) k X(o) ~ g7(o)«i^ ) (o)fc- D (o)i( & - D (o) fc X(o)^(o)'x(o)) 
-^ a (o)(bdiX(o)djX(o)) +-D(o)(i^(o) fc (^ fe X(o)^)X(o)) - ^(o)i^(o)#^(o)fcX(o)^(o) fe X(o)) 
+ g7(o)*i a (o)(^(o)fcX(o)^(o) fc X(o)) + & - D (o)i- D (o)iX(o) n (o)X(o) ~ - D (o)(i(W(o)X(o)- D (o)j)X(o)) 
+^7(o)ii^>(o)fc(^(o)X(o)^>(o) fc X(o)) - |7(o)ijK D (o)X(o)) 2 



+c ^27(o)ii £) (o)fcX(o) £, (o) X(o) - 2 ^X(o)^'X(o) J D (o)iX(o) D (o) X(o) 
+4k 2 ^7(o)ik7(o)ji - ^7(o)fcz7(o)ij) — :p=%) fe '- (3-19) 



25 



Similarly, the first order equations for the coefficients in the dilaton expansion are 

A(o)-r 2 A(o)^(A( )) = 0, (3.20) 



A( 2 ) - + 2r 2 A (0) t/ / (A (0) ) + r 2 A(o) 2 £/"(A (0 ))) A (2) = C 2 \o) 2 Z\\ ))Rb( )] 
-M(A ( o))D(o)A ( o) + ^(q)- 1 (2M(A (0) ) - A (0) M'(A (0) )) D (0)i X {0) D 
+^r 2 A(o) 2 e , (A (0 ))Z)(o) l X(o) J D(o) J X(o), (3.21) 



, 2 

\2) 



A ( 4) - + 2r 2 A ( o)^' + r 2 A( ) 2 C/") A (4) = (u' + 2A (0 )C/" + ^A (0) V") T 2 A ( 
-A( 2 )M / D( )A(o) + M ^7( 2 )^L>( )j J D(o)iA(o) - □(o)A(2) + 2( 2jD (o)j7(2)j - D (o)n{2)) D {o) j \$) 

-2 

2~A( )A( 2 ) (2©' + A( O )0") -D(o)jX(o)^(o) l X(o) + A( 0) A (2 ) (2S' + A( )S") i?[7( )] 



r 2 



A(o) 2 @' (2 J D(o) i X(o)-D(o)iX(2) - 7(2) y ftx(o)5iX(o)) 



2 

+C _2 A(o) 2 " (D^D^j^ij ~ D (o)7(2) - 7(2) y -R[7(o)]ii) 

+^A (2) A (0) ~ 1 (M' - A (0) M") £> (0)i A ( o)iVA(o) + ^ A (o) _1 ( M ' " 2A (0) ' 1 M) x 

(A(2)- D (o)jA( )-D(o) J A(o) + A(o)7(2) y ^A(o)9jA( ) - 2A(o)-D(o)iA( ) J D(o) l A( 2 )) 

i' (%[7(o)]^[7(o)] " ^ 2 [7 ( o)]) +c' (D { o)a(o)^(o)^(o)) 2 
^•[7(0)] - ^#[7(0)b(o)«) £(o)*X(o)£(o) J X(o) - ^( D (o)X(o)) 



+ 4^2 A (0) 



+6' 



+ K 2 r 2 A( ) 2 -^^A(4)- (3.22) 
Finally, for the axion we get 

X(0) = 0, (3.23) 



X( 2 ) - 7X(2) = -Z~\\ m ) (e(A ( o))D(o)X(o) + ®'(\o))D {0 )i\o)D { o)%o)) , (3-24) 



x( 4) - |X( 4 ) = -A( 2 )^ 1 (e' - z~ l z'e) u mm + z^e^'w^D^x^ 

-Z~ 1 ®' (-D(o)iX(o)^>(o) l A(2) +^(o)iX(2)^(o) l A( ) - l(2) lJ di\o)djX(o)) 

-z- l \ {2) (e" - z-^'e') D^Aioi^toi'xfo) 



2 Z D (0) l \b yRij[y(p)] - g-R[7(o)]7(o)«J ^(0)^(0) + ^(oj^^fojXfo)) 

+2cD( )fcX(o)^(o) fe X(o)^X(o) ) +k 2 ^ -1 — — 5r X (4)- (3-25) 

/ v7(o) 
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In all these first order equations all functions are functions of Am) (r, x) and primes denote derivative 
w.r.t. A( ). 

From these equations we can now easily construct the asymptotic expansions. The non-linear 
equations for the zero order coefficients can be integrated exactly in terms of the function U(X^). 
Using the asymptotic form of U (A) in (|3.4|) one obtains the asymptotic expansions 

ATT \ 

+ 3(a-2)6g [rj + ° (WW*)/*" 2 ) 1 9m&> ( 3 - 26 ) 
+0 (log 2 (r/£)/r 3 ) , 

X(o)( r ^) = X(o)0»0> (3-27) 

where g(p)ij(x), X^(x) and X(o)( :c ) are arbitrary and we identify them with the sources of the dual 
operators. A striking feature of the zero order dilaton expansion is that the source, A( )(x), that 
couples to the dual operator appears not in the leading order, but at 0(l/r 2 ). This is related to 
the fact that \rg)(x) can be removed by a bulk diffeomorphism corresponding to shifts of the radial 
coordinate, which we will demonstrate below. However, we have already made extensive use of the 
identity d{\ = 0(X 2 ) in deriving the results in Tables [5] and [H which follows precisely from the 
observation that the source of the dilaton appears at subleading order. 
At the second order we get 

l{2)ij{r,x) =g(2)ij{x) + f~log(r/i) 

/A 2+252/3 (- /I 
A( 2 )(r,x) = (^-J (A (2) (*) + <^-log(r/^ 

X(2)(r,x)= (;)^ /3 (x(2)(^) + o(ilog(r/^y (3.28) 



where 



I 2 ( 1_ 

9(2)ij(x) = -— ( Rij[g(o)} - ^9(o)ijR[9(o)} (3-29) 
- ( M pi N c) 1 Zo (diX(o)djX(o) ~ ^9(o)ijD(o)kX(o)D(o) k X(o)Y) , 

A(2)0»0 = -j (J^ R \9(o)] -n ( o)A (0 ) (3.30) 

+ ^2 (^1 " ^(0)^(0)^(0)^(0)) , 



r/2 



X m (x) = ^(M 3 iV 2 ) 1 D (0 )X(o)- (3.31) 
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Finally, at fourth order the asymptotic expansions take the form 

/£\ 2£ 2 /3 // r x_ /log(r/f) 
7(4)ij(r,x) = ( - J 9(4)ij(x) + log(r/£)g (4)lj (x) + g {4)lj (x) + O ( 

A(4) (r, x) = (;) ^ ^ (a (4) (x) + O (^tl^J ) , (3.32) 

X(4)(r,x) = /3 ((-) X(4)(x) + log(r/£)x (4) (x) + x w (x) + O (M^l 

where the terms gu)ij(x), A( 4 )(x) and X(a)( x ) are undetermined and are therefore identified with 
the normalizable modes. The flow equations relate these to the undetermined renormalized mo- 
menta, thus leading to the expressions for the one-point functions in terms of the coefficients of the 
asymptotic expansions, which we present below. 

In order to write down the explicit expressions for the coefficients in these fourth order expan- 
sions, it is useful to define the traceless tensors 



1- \ ^kl 2„/„ 1 



Huj — —2 yRikji — -^Rki9(o)ij J R + g-R ( Rij — -^Rg(o)ij 

+-^D^iD^jR - O^Rij + ~<7( )j,- □(())#, (3.33) 

H 2 ij = (—Rikjl + -j9{0)ijRki - R(ik9(0)j)l 

1 1 1 \ 

-QR9(o)ij9{o)ki + -^Rg{Q)ik9(Q)ji + ■^Rijg(o)ki) X(o)-°(o) X(o) 



+D(o) {i D (0) (d k X(o)dj)X(o)) ~ 2 D (o)(^X(o)<9iX(o)) 
-Qg(o)ijD(o)kD( )i{D(o) k X(o)D(o) l X(o)) 

-g-D(o)i-C»(o)j(-D(o)fcX(o)-C'(o)' c X(o)) + Qg{o)ij^(o)(D{o)kX(o)D( ) k X(o)) 

+ D (0)i D (0)jX(0)^(0)X(0) ~ D (0)(i(^(0)X(0) D (0)j)X{0)) 

+ Qg(o)tjD(o)k(0 (0) X(o)D {0) k X(o)) ~ ^5(o)ij(n(o)X(o)) 2 , (3-34) 

H 3ij = (^9(o)ijD(o)kX(o)D(o) k X(o) ~ 2 diX(o)djX(o)^J ^(o)*X(o)^(o/X(o)- (3.35) 

As mentioned in the previous section, these correspond respectively to the derivative of the three 
terms proportional to ci, c 2 and C3 in the boundary term (|3.8p w.r.t. the induced metric. All 
curvatures here are curvatures of the boundary metric gm)ij- The first two coefficients in the fourth 
order expansion of the metric are just linear combinations of these three traceless tensors. Namely, 

g m (x) = jH Uj - jZ {M^)- 1 H 2ij + ^Zl (Af* N 2 C )~ 2 H m (3.36) 



g{A)ij{x) 



t_ ( 2U 2 
12 \bl 



lij 



4 Z ° 



nu 2 2g 2 z 1 \ 

V3&o 3 b o z o) 



2i; 



+-Z 2 (M^N?) 



-2 (m 

V36g 



2£ 2 



+ 



2Z 1 
boZ 



H- 



3ij- 



(3.37) 
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Moreover, for the coefficients of the axion expansion we have 



'(0) 



my 1 

r3 Ar2\-2 



R ij[9(o)] ~ o R \S{o)]9(o)ij) -D(o) J X(o) + 5-D(o)i n (o)X(o) 



-Z (M^N*) D {0)k x {0 )D {0) k X(o)diX(o) , 



(3.38) 



and 



X(4)( 



-D 



(0) 



3 + 6 -^o 



2 {8U 2 2£ 2 2Z X \ - n ft- , 



3 + b Z 



(3.39) 



Finally, we can write down the general expressions for the one-point functions, i.e. the renor- 
malized momenta, in terms of the coefficients of the asymptotic expansions. Starting from the 
dilaton, the exact one-point function in given by 



{Ox) 



^2" { — £^ A (4) +Rij[9(o)\R J [9(o)\ ~ 3^ [9(o)\ 



-2Z (MfrN^ 



1 1 

R ij[9(o)] - g-%(o)]0(o)ij 1 ^>(o) J X(o)-0(o) J X(o) - 2( D (o)X(o))' 



+^Z 2 (M p 3 ,JV c 2 ) 2 (5(0)^(0)5(0)^(0); 



(3.40) 



Notice that this gets contributions from the dilaton normalizable mode, A(4), plus a combination 
of the metric and axion sources which is nothing but the conformal anomaly of the axion-gravity 
system in strictly asymptotically AdS space (cf. (|B.12[) ). 

All results quoted so far are valid for arbitrary dilaton source A( )(x). However, in order to 
simplify the expressions for the stress tensor and axion one-point functions, we will only give the 
expressions for constant \(p), independent of the transverse coordinates. This is not a big disadvan- 
tage since, as we shall see, the dilaton source can be removed or restored by a bulk diffeomorphism 
corresponding to a Weyl transformation of the boundary metric. The full dependence on a generic 
dilaton source can therefore be restored starting from the expressions given below for constant 
A(o)(x), by a suitable boundary Weyl transformation. 

The one-point function of the stress tensor can be written in the form 



(Ti 



l l I r-en 



2 1 

— (fiy - TcQ,g m ) - — {O x ) ren g {0)ij , 



(3.41) 
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where the tensor fijj is given by 



n 



ij - 9(4)ij + — g-A( 4 )fli(o)ij 



-j^Huj + jZ {M^N 2 ) 1 c' 2 H 2ij - ^Z 2 (M^Nt) 2 c^- 

^2 / _ ]_ _ j 

+-J ( D(Q)kD(o)(i9(2)j\ — -^((WQoma — —D(n\iDin\ iTr i 



)j) - 2 U (0)5 , (2)ii - 2^(0)i^(0)jJ-r 5(2 ) 



~9m ( R + Z o{K N cY l D ( o)kX(o)D (0 ) k X(o)) 

I 2 ( - - — 
~^9(o)ij {D( ) k D {0) l g-( 2 )ki ~ □(o)Tt3( 2 ) - 9(2) kl R, 

I 2 -if- 1 

-— Z (M^N 2 ) g (0)lj ( D (0) k X(o)D (0)k X(2) ~ ^9(2) H d k X(o)di 

t 
192 



7X(0) 



Rki - ^R9(o)ki D {0) X(o)D( ) X(o) - ~( D (o)X(o))" 



-2Z {M^N 2 ) 



+-Z 2 (M^N 2 y (5 (0) l X(o)5 (0 )iX(o)) ■ 



1 — ^ - fc-. f> i - 1 ^— - ^2 



Here we have defined the shifted constants 



2U 2 (2£ 



i (\e 



c 'i = c i + + Hr " 1 " 2 Clog6 + 7 Hr ~ 1 



3^ 

4 = C2 + ^+(^-l)6 A (0) -(2C 

C 3 = C3 + ^+f^- i ).V,.„ ^ (2C 



2Zi 



1 /4£ 2 



1 /4£ 2 
Iog6o + T 



(3.42) 



(3.43) 



In principle, one can set these constants to zero by suitable choice of scheme, i.e. by a suitable 
choice of ci, c 2 and C3, but we give the full expressions so that one knows exactly what scheme 
needs to be chosen to achieve this. 

Finally, the one-point function of the operator dual to the axion is 



in \ 4Z ° ~ 

\ U X/ren ~ _ T27 X ( 4 ) 



K 2 £ 



k- 



-^Z 2 (M^N 2 ) 4D (0) * [D m x {0) D (0) k X(o)d t X( ) 



3 



+ 2^ Z ° 

n (0 )X(2) - « (25 (0 )i^(2)5 - £>(o);Tr5 (2) ) D (0 f X (o) ~ 9{2) 3 ' D (0 )iD( )jX(o) 



(3.44) 
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where again we have introduced the constants 



° 2 = C2 + 36§ + — 6oA W " W 



2U 2 , 2^ 
3 



C 3 = c 3 + ^if + — & oA( ) - 



36g 



(3.45) 



to abbreviate the above expression. 

3.3 Asymptotic diffeomorphisms and Ward identities 

Now that we have determined the general form of the asymptotic expansions for IHQCD and we 
have identified the exact one-point functions, we can proceed with the derivation of the holographic 
Ward identities. These follow as a consequence of the existence of a class of asymptotic bulk 
diffeomorphisms that preserve the structure of the asymptotic expansions. 

Let us consider a generic infinitesimal bulk diffeomorphism, 5x^ = — and demand that it 
preserves the gauge fixed form of the metric, namely that it does not modify the lapse and shift 
functions. This requirement leads to a pair of equations for the vector field namely 

-CgSrr = C = 0, 

C,H,, -,j^ J ■ <> J C) <>• (3.46) 
where £g is the Lie derivative w.r.t. the bulk vector Solving these conditions gives 

C = Sa(x), 

poo 

C = + djda(x) / drW, x), (3.47) 

J r 

where u{x) is an arbitrary function of the transverse coordinates and Co( x ) is an arbitrary transverse 
vector field. Inserting now the asymptotic form of the induced metric we obtain 



e = &*) + ^ (f 1 



-2r//_ -2|— 1 



log(r/£) 



(3.48) 



Under this bulk diffeomorphism then the induced fields transform as 

2 



S aij = £ i9ij = L aij + 2KijC = H olii - —UiX^jdaix) + / 3 ), 
5^X = £ 5 A = L^X + eX = ejfltA + r 2 A 2 ^W) + 0{r-^l\-^l% 
kx = Ct-X = hx + ex = CodiX + 0(r- ' 



~-2?73 e -2r/ 



(3.49) 



where denotes the Lie derivative w.r.t. the transverse components of the bulk vector field £. 
It follows that the sources, g(p)ij(x), A( )(x) and X(o)( x ) transform under such a diffeomorphism as 



%(0)ij = 


L Zo9(0)ij + 


j5a(x)g {0)ij , 


k\o) = £ 


l o (x)di\ 0) - 




^X(o) = * 


Z(x)diX(o)- 





(3.50) 
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There is nothing surprising about the transformation of the metric and axion sources. They are 
exactly as they would be in the case of strictly asymptotically AdS space. Namely, the bulk diffeo- 
morphisms that preserve the form of the asymptotic expansions correspond to arbitrary boundary 
diffeomorphisms parameterized my d(x), as well as boundary Weyl transformations, parameter- 
ized by the function a(x). What is rather unusual, is the transformation of the dilaton source, 
A(o) under the Weyl transformation 5a. Contrary to the usual multiplicative transformation of the 
sources, the transformation of A( ) is additive under boundary Weyl rescalings. This means that 
one can in fact remove the dilaton source completely by means of a boundary Weyl rescaling. 

We can now also determine the transformation of the one-point functions under boundary 
Weyl rescalings. The above transformations of the sources imply that the corresponding functional 
derivatives transform as 



J9(0)ij 

Moreover, the renormalized action 



-5a(x)-^- — , 5 a 



i.) 



SX 



(0) 



o, s a 



0. 



lim S( 

r—too 



(3.51) 



(3.52) 



is invariant under any bulk diffeomorphism since the boundary term (I3.8P does not break the bulk 
diffeomorphisms. Hence, 



-M( T i) > 

\ J I ren 1 




5a (O 



XI ren ' 



(3.53) 



that is all one-point functions transform homogeneously under boundary Weyl rescalings. 
The Ward identities now follow from the identity 



1 



d d x -^g (0)ij (ri) ren + 5fX {0) (O x ) ren + 8&(o) (O x ), 



0. 



(3.54) 



Inserting the transformation of the sources under the bulk diffeomorphisms considered above and 
using the fact that 5a(x) and Co( x ) are arbitrary leads respectively to 



" * / ren 



^ i^^ren ' 



D W (Tjlen + ^ren fy\o) + ^ X ) ren ^X(O) = 0. 



(3.55) 



An immediate consequence of the trace Ward identity is that the tensor introduced above is 
traceless. 

Finally, let us examine a bit closer the relation between the dilaton source A(q) and boundary 
Weyl transformations. Note that the assignment of sources in the leading order asymptotic solutions 
(I3.26P is rather arbitrary. In particular, we could have included a factor of the dilaton source in 
the definition of the boundary metric. Suppose, in particular, that we define 



9(o)ij(x) 



(3.56) 
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Then, inserting this boundary metric back in and evaluating the variation of the induced 

fields with respect to variations of the dilaton source Am) , one immediately sees from (|3.49p that 
transformations of the dilaton source correspond precisely to boundary Weyl transformations upon 
the identification 

* ' (3.57) 



a{x) 



-MA( )( 



This observation proves that the source of the dilaton is gauge freedom that can be removed by a 
bulk diffeomorphism corresponding to a boundary Weyl transformation. This property can also be 
used to restore the full dependence on the dilaton source of the one-point functions of the stress 
tensor and the axion. 



4 Concluding remarks 

We considered a generic dilaton-axion system coupled to Einstein-Hilbert gravity in arbitrary space- 
time dimension and we carried out the procedure of holographic renormalization of this action for 
dimension up to and including five dimensions. The general boundary term that renders the varia- 
tional problem for this action well defined is summarized in Tables [2] and [3l This result is applicable 
to a very wide range of holographic models in the literature, including N = 4 super Yang-Mills in 
four dimensions, Improved Holographic QCD and non-conformal branes. We explicitly evaluated 
this general boundary term for a constant dilaton potential, corresponding to the standard dilaton- 
axion system dual to the complexified coupling of N = 4 super Yang-Mills in four dimensions, in 
Appendix [Bj and for IHQCD in Section [3j In particular, we systematically derived the generalized 
Feffer man- Graham asymptotic expansions, provided exact expressions for the one-point functions 
in the presence of sources, and proved the holographic holographic Ward identities by studying the 
asymptotic bulk diffeomorphisms that preserve the form of the asymptotic expansions. 

In the case of IHQCD, an important lesson from the analysis is that the source of the dilaton is 
not a physical coupling, but its value can be thought of as an energy scale. In particular, changes 
in the dilaton source can be absorbed by a Weyl rescaling of the boundary metric. Moreover, the 
operator dual to the dilaton field A is the combination 

0A = /?(Ayjv/)TrF 2 , (4.1) 

which has fixed scaling dimension 4 under renormalization group flow. This is the combination 
that appears in the trace Ward identity 13. 55^ since the coefficient relating the trace of the stress 
tensor to the operator 0\ is a constant, independent of the renormalization group scale A( ). 

Our calculation for IHQCD, independently of the legitimacy of the model as a physically sound 
holographic model dual to pure Yang-Mills theory in four dimensions, provides us with an explicit 
example of a gravity model that can accommodate operators with running scaling dimensions. This 
is particularly interesting since it allows us, in principle, to develop supergravity holographic models 
capturing the dynamics of operators with non-protected scaling dimensions. 
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Appendix 



A Functional integration 

In this appendix we outline the derivation of the functional integration formulas in Table [TJ In 
particular, the question we want to address is the following: given a local functional TZ(2n){ ( f) of 
the scalar field, (p, what is the local functional F(2 n ){^P) such that (cf. (|2.38p ) 

j^e a ^K (2n) (cp) = VW^) + e ani ' p) d i v (2n) i (ip,5ip), (A.l) 

where a n ((p) is a prescribed function of tp and t>( 2n ) l ((/?, <5</?) is some vector field? If the source 
7^(2n)(v) does not involve spacetime derivatives of the scalar field the answer to this question 
is given by simple integration. However, if TZ^njif) involves derivatives of <p, then determining 
F(2n){ ( p) becomes less trivial. It turns out that one can still find a general formula if 7£(2ra)(v) 
involves first derivatives of the scalar field, but once Tl(2n)( i p) contains second and higher order 
derivatives of the scalar field finding a general formula becomes much harder. What we will do 
instead here is to consider only the sources Tt(2n){ l f) which are relevant to our computation of the 
solution of the Hamilton-Jacobi equation in the main body of the paper. 

• ^(2«)(v) =r 1 m((p)t ili2 ~ im d il <pd i2 <p...d im <p 

The first example is a generic source that is polynomial in first derivatives. Here, t lll2 '" lm is 
an arbitrary totally symmetric tensor that does depend on ip. In this case we can write F(2n) 
as 

F { 2n) = e a " (a(<p)t^-- A ™d ilV d i2 <p. . . d im y + D { {f3^)t ti ^^d l2 ip . . . d im y)) , (A.2) 

where the functions ct((p) and /3((p) are to be determined. Evaluating the variation of this 
expression and inserting the result in (1A.1|) we obtain the two equations 

a m 
p = —a 

a' 



■ ™ a 'n . r l 



c/-«^-fr, (A.3) 

which can be solved to determine a(ip) and j3((p). Since /3(</?) contributes a total derivative 
to Hamilton's principal function, i.e. to e~ a " n ( tp > F( 2n ), we are only interested in a(<p), which 
is given by 

aft*) = aTe~ a " [* ^ e ^ m r lm (0) = f ri m(<p). (A.4) 



n.rn 



Similarly, for a source with a single second derivative of the scalar field we can write 

F { 2n) = e an (a(<p)t ij DiDjtp + p(<p)t ij dupdjtp + A [j(v)t ij d jip + 6(<p)Djt ij )) . (A.5) 

Here, t %3 is again a symmetric tensor that does not depend on <p. Inserting the variation of 
this expression in (|A.1|) leads to the following equations for the functions ot(tp), /3((p), 7(9?) 
and 5(ip) 

a'Jfi + 7') + a" - = 0, 
a' n ( 1 + 5')+2(a'-P)=0, 
a' n S + a = 0, 

2a / + aU« + 7)-2/3 = ^. (A.6) 
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These can be immediately solved to obtain 

a(<p) = a' n e- a - f ^^rtf) = f r{(p), (A.7) 



n,l 



U' 

/%) = -die— f ^a'- 2 U'a' n dl (£) a {0) = - £ U> J n d% (i-) a(0). 

(r 1 2 2 (<p)t t ^ kl + s 1 2 2 (y)*2 W ) di<pdjipD k Di(p + (r 2 2 (^)t^' W + s 2 2((p)t l 2 jkl j DiDj(pD k Di(p 

As a final example we consider a generic term with four derivatives, but we restrict to co- 
variantly constant tensors t^ kl and t% 1 ■ In particular, for our purposes it suffices to take 
these two tensors to be the two linearly independent tensors constructed out of the metric. 
Namely, we will take 

qkl = 1 ^ifcyj + yiyfc + ^ qkl = 1 ^ifcyi + yjyfc _ g^J^ . ( A .8) 

Moreover, we need not consider a source for four first derivatives since we have already 
computed the result for an arbitrary number of first derivatives above. Writing then 



F(2n) = e a " [AWDiDjvDkDvp + B^ kl d iV d jV D k D lV + C i i kl d i <pd j <pd kV dnp 

+ A (yE^ kl DjD k Dnp + H^ kl d jV D k D lV + G^d^d^d^ ) , (A.9) 
and inserting the variation of this expression in (|A.1|) we obtain the set of coupled equations 
3A njki + a ^ A ijkl + B iikj + B ikjl _ 2B ijki + a > nII ijkl = _L (r 22 ( (p ) t » kl + s 2 2{ip)-ei kl 



2 A >ijki _j_ 2 ( B >kju _j_ gtljik \ _ -yic 1 ^ 1 + a' ( B^ kl + H H i kl + G kli ^ + G l -* ik + G k i 



id 



— (r 1 2 2 (<f)t\ 3kl + s 122 ((p)t^ kl 

_ 3(J Hjkl + B n(ijkl) + ^ (fjijkl + G /(ijfcO^) = 0, 

^ A 'ijki _|_ ^fciji _|_ ^?jfc _ 2B^ kl + a' (^E'^ kl + H^ kl " 
2A ijkl + a' n E ijkl = 0. (A. 10) 

Here, the parentheses in the indices mean total symmetrization. Note that these are five 
equations for six undetermined tensors. They can be solved as follows. First we can use the 
last two equations to eliminate E l i kl and JEP 3 . This leads to the decoupled equation 

A ' ijkl + [ a ' n ~ <9 Aijkl = h fa^" + **(*o# w ) , (A.n) 

whose solution is 

A*» = c4V- f ^a>- 2 (r 22 (<p)tf l + s^^tf 1 ) . (A.12) 

Next, we can set 

a > n Qijkl = 4C ijkl^ ( A13 ) 

so that we obtain the following two equations for B^ kl and C^ kl : 

2 B >ijkl _|_ B 'kijl _|_ B llijk _|_ Z^/ _ 2°n \ ^yfci _j_ 2a n / ^fcjj'Z _|_ B lijk\ 

\ n < ) < K - ' 

= ^ (ri> 2 (^ W + W^f ') - K<? (£) (A.14) 
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and 



C Hjkl + (a' n - ^pj C ijkl + B"^ = 0. (A. 15) 

The last equation for C^ kl gives 

C** 1 = a'*e- a " J* ^e a "a'~ 4 (-B"& kl ^ . (A.16) 

To solve for B^ kl we first decompose it as 

B^ kl = a(ip)t\ jkl + (3(<p)t 2 jkl . (A. 17) 

Inserting this into the above equation for B^ kl leads to the two decoupled equations 

4a' + a' n a = —a, 



where 



co = sMv) - 2U'a' n % (J^J a'*e- a " ^e a -a'~ 2 s 2 ,(Cp). (A.19) 



These equations give 



P = f^e a "a>- 3 u. (A.20) 

Moreover, since B^ kl ^ = at^ kl , we have 

Cijkl = _ t W a >4 e -an ^ e a„ a /-4 a // (A 21) 

These results can be simplified considerably by noticing that the terms involving a combine 
into a total derivative, up to a homogeneous term that is irrelevant since it contributes a 
finite piece. To see this, integrate by parts the first term in 

at^ kl d i ipdjipD k D l ip + ct^ kl diipdjipd k ipdiip, (A. 22) 

which gives 

— 2atj' kl diipdj(pDkDiip + (c — a')t^ kl diipdjipd k ipdiip. (A. 23) 

We can now replace the first of these expressions in Fr 2v \ by 2/3 of the first expression plus 

1/3 of the second, such that the coefficient of t l f ki 'diipdjipD k Diip vanishes. This replaces c by 
c = c — a! /3. However, we have seen above that the only source of the equation satisfied by 
c, or c now, is the coefficient of t l ^ kl diipdjipD k Dnp. Since we have now set this coefficient to 
zero, c satisfies a homogeneous equation and hence we can set c = 0. We therefore conclude 
that, without loss of generality, we have 

B^ kl = (3(ip)t 2 jkl , C ijkl = 0, (A.24) 

for any a (if). Notice, in particular, that the source r 1 2 2 (ip) gives no contribution at all. 
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B Dilaton-axion system with constant dilaton potential 



In this appendix 13 ! we focus on the special case of a constant dilaton potential (£ = 1) 

V(tp) = d{d - 1) = 12. (B.l) 

This special case drastically simplifies the solution of both the zero order problem (|2.26p and the 
linear recursion equations (|2.3ip . In particular, the first order differential equations are replaced 
by algebraic equations. Namely, we havj^l 

U 2 + d{d - 1) = 0, (B.3) 



d- 1 
d-2n 
d-l 



U£ (2n) =K {2n) , n>0. (B.4) 



From these we immediately obtain 

U = -(d-l), 

£(2) = - 2{d l 2)K 2 Vl {Rh} - di<p&<p - Z{ V )d lX d l x) , (B.5) 

where we have fixed the sign of U by demanding that the solution leads, via the relation 

2 

to the correct asymptotics for the induced metric jij. 

In order to compute £(4) we need to compute the momenta from £( 2 ) • We easily get 

= 2{d-2)n^ ~ " Z ^) dt X&X ~ \l 13 (R ~ d k cpd k <p - Z(v)d kX d k X 

*x&) = - {d ^ 2 )K^ ( Z ^ n ~> x + Z'(<P)ai<P&x) ■ (B.7) 



13 The author is grateful to David Mateos and Diego Trancanelli for checking the results of this appendix and for 
pointing out typos in a preliminary version. Of course, the author is solely responsible for any remaining typos. 

14 In addition to the constant solution U(tp) = — (d — 1) of (|2.29|) there is additionally a one-parameter family of 
non-constant solutions given by [301 131] 



' ( r -l - -(,/- li<-sh j ^j-J-^-ipo) J ^ (B.2) 

for some arbitrary constant <fi and AdS asymptotics requires that tp — > ip asymptotically. However, this solution 
only allows for a constant non-normalizable mode for the dilaton, i.e. ip a cannot be a function of the transverse 
coordinates x and hence it does not correspond to the most general asymptotics. In fact, a domain wall with such 
a 'fake superpotential' describes a vacuum where the operator dual to the dilaton field has acquired a VEV |31j . 
Finally note that a special feature of equation (I2.26|l with constant scalar potential is that the solution describing the 
most general asymptotics, i.e. U(tp) = — (d— 1), is isolated in the space of solutions, while the one-parameter family 
of solutions describes special asymptotics. This is the reverse of what happens for non-constant scalar potentials. 
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Hence, from (|2.32p . we evaluate 

= -^2(^)2^7 (R ij Rij + (d^d^) 2 + z 2 ^)(d lX d l x ? - 2&a i <pd j < P 

-2Z{v)R ij d iX d jX + 2Z(<p)(d i <pd i <p)(d jX d j x) - ^^( R - d ^ d ^ - Z (v)d lX d l x) 2 
+ (n j tp-±Z'(v)d iX d i X ^ + Z^)(D jX + d tp logZ(ip)d^d i X ) 2 j. (B.8) 

Now, the recursion relations tell us that £( 4 ) is given by 

(d-4)£ (4) =K {A)1 (B.9) 

which is ill defined in d = 4. This problem is of course well known and it is related to the breakdown 
of full diffeomorphism invariance of the Hamilton-Jacobi functional, which in turn leads to the 
conformal anomaly of the dual conformal field theory [T5] . It was shown in [T7J that the Hamilton- 
Jacobi approach can be still applied in this case and reproduces the results of the Fefferman-Graham 
expansion [B] provided the radial cut-off is related to the deviation of d from the desired value of 
4. In particular, we set 

r„ = (B.10) 

and define £(4)|r = — 2^o^(4)|r > so that 

£ (4) =-^ (4) . (B.ll) 

This then leads to the full counterterm action for d = 4. Namely, dropping again the subscript 
from the radial cut-off, 



Set = d A x^j ^3 + i (R - dupd'v - Z{!p)diX&x) 



+2Z(<p)(d i <pd i <p)(d jX &i X ) - ±(R - dw&tp - ZMdiX&x) 2 

+ (n y <p - ^Z'i^dixd'x^j + Z(<p) (D 7X + d v log ZMditp&xfj J • (B.12) 

This expression does not quite agree with the one given in (26) of [32], which we believe is incorrect. 
B.l Asymptotic expansions 

The above recursive solution of the Hamilton-Jacobi equation, implies that the canonical momenta 
take the form 

7T* J ' = 7T {0 f + 7T (2 f + (-2r)7f (4 ) l ' j + 7T (4 f + ..., (B.13) 

and similarly for the dilaton and the axion, obtained by differentiating Hamilton's principal function 
w.r.t. the corresponding induced field. The form of these expansions implies in turn that the 
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induced fields admit the asymptotic expansions 



e (S(o)ij + e g {2 )ij + e 



2r/t( 4 )jj +5(4)ij) + • • •) 
2r^ (4) + V9 (4) ) + . . . , 

X = X(o) + e X(2) + e (- 2r X(4) + X(4) ) + ■■■, 



V = V(o) + e 2 V(2) + e 4r 



(B.14) 



where the dots denote subleading terms that are unambiguously determined in terms of the terms 
shown. These are precisely the well known Fefferman-Graham expansions 

The coefficients in these expansions can be easily deduced from the expressions for the canonical 
momenta we obtained in the previous section. In particular, inserting the asymptotic expansions 
for the induced fields in the expressions (|2.6p for the canonical momenta on the one hand, and in 
the above expansions of the momenta on the other and comparing the two determines at second 
order 



<P(2) 
X(2) 
9(2)ij 



+ 



1 



□(o)^(o) - 2 Z '(^(o))2(o) y diX(o)<9jX(o) 
Z'(ipto)) 

D (o)X(o) + ^— j 9(ofdiX(o)d3<P(p) 

\ ( R [9(0))ij ~ di¥>(o)<9j¥>(Q) - Z(<P(o))diX(o)djX(o)) 
^9(a)ij (R[9(o)} ~ <9fc¥>(o)dV(o) ~ z {V(o))dkX(a)d k X(o) 



(B.15) 



Comparing the logarithmic terms gives 



1 



1 



K 



lim 



:7T t/ 



r—too ,y^y f (4) 

4 „, ,. 1 



r— s>oo 

hm — — 

r->oo ^y Slfi 



lim 



1 



— {h{ij j -Trh (i) g (0 f) = lim — 7f (4) J - 



(TxC 



(4), 



<FxC 



(4), 



— Z(^ (0) )x (4) = r lhn — 7f v(4) 



r 1 6 
hm — — 



dxC 



(4), 



(B.16) 



which can be easily evaluated explicitly using the above expression for but we will not need 
these expressions here. Finally, the 0(e~ 4r ) terms lead to the expressions for the renormalized 
one-point functions in the presence of sources, namely 



(O 



ip/ren 



-7T, 



mo) 



V>(4) 



4 2 1 

"72^(4) - 72^(4) + ^2 (-5(2)* J ^(0)i^(0)iV(0) 



1 



D( ) t 9(2)ij ~ 2 L> (0)i Tr f(2)J -D(0)V(0) + D (0)^(2) 
■^"(^(0))^(2)5'(0) iJ ^X(0)^X(0) + ^'(^(o))5(2) y ^X(o)5iX(o) 



-Z'(^(o))ff(o) iJ <9;X(o)<9jX(2)) , 



(B.17) 



X/ren 



-7T 



4 2 2 

Uh = -^2 Z (^(0))X(4) - ^O(0))X(4) - ^^(V(0))V(2)X(2) 

+ ^~2 ( Z '(^(0))V(2)n(0)X(0) - ^(^(0))5(2) U ^(0)i^(0)jX(0) 



D(p)g$)ij ~ 2 D (0)J Tr 9(2) J D (o) 3 X(o) + Z(¥J(o))D( )X(2) 

+^"(V(o))^(2)'9iX(o)'9V(o) - £ / (P(o))fl , (2) y fii¥>(o)fyX(o) 
+^'(>(0))diX(2)<9V(0) + ^ / ( ( ^(o))^X(o)^V(2)) , (B.18) 
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2 2 1 

(Tij)ren = -=^{A)ij = ~2 (fl(4)ij ~ ^9(4)9(0)^) + ~2 (^(4)ij ~ Tr /i( 4 )flr( )t.j) 



1 



"4 (^5(2)5(2)^' - Tr (5( 2 ) 2 )5(o)ij) 



+ 2^ [ D (0) kD <PKi9<$3) ~ 2 n (0) 5 ( 2 )^' ~ ^W^Wi 1 ^ 2 ) 
-2%¥>(o)<9,-)<£( 2 ) - ^'(y(o))v 3 (2)5iX(o)5jX(o) - 2^(v2(o))%X(o)^i)X(2) 



5(2)ij (^b(o)] - d k ¥(o)d kl P(o) ~ z {<P(o))dkX(o)d h X(o) 

1 



+^9{o)ij9(2) kl (-%(o)]fcz - d k ip(o)diip(p) - Z(ip {0) )d k X(o)diX(o)) 
;5(o)« ( D (0)kD {0) ig(2) H ~ 0(0) ^ 5(2) - 2<9 fc ^(o)dV(2) 



2' 
1 

~2~' 

-^(^( ))^(2)^X(0)S fe X(0) - 2Z(^ (0) )a fc X(o)^X(2))) • (B.19) 

Note that in some places we have symmetrized indices with weight one, i.e. (ij) = ^(ij It 
can be easily checked that the expression for the one-point function of the stress tensor evaluated 
at zero dilaton and axion sources completely agrees with the expression in (3.12) of [33jP^I 

The holographic Ward identities can now be derived as in the main body of the text in the case 
of IHQCD. Namely, the second of the equations in (I2.18P holds order by order in the expansion of 
the Hamilton-Jacobi functional in eigenfunctions of the dilatation operator. In particular, we have 

D (0){Tij) ren + (O v ) ren D(p)jip(ti) + {O x )renD(fi)jX(0) = 0. (B.20) 

Moreover, considering the transformation of the renormalized action under the asymptotic diffeo- 
morphisms that leave the Fefferman-Graham expansions form-invariant leads to the trace Ward 
identity 

(T-} ren = A(g(Q),ip( Q ),X(p)), (B.21) 



where 



,4(-. T r. x) ^ (wRi, + + Z-\<p){d iX d l xY - Wdwdjtp - 2Z(<p)RVd iX d jX 

d 



in 2 \ 



+2Z(<p){d i <p&<p){d jX &'x) ~ j^—^ ( R ~ ~ Z^dad'x? 

+ (u^ip - ^Z'{ip)diX&x) + Z{<p) (D 7 x + d v log Z{<p)di<p$x) 2 \ , (B.22) 

is the conformal anomaly. 

As a final comment, we observe that since 

£ (4) = = ±yfyA, (B.23) 

the terms involving huuj, (p^ and X(4) i n the above expressions for the renormalized one-point 
functions can be ignored since they are scheme dependent and they can be removed by adding the 
finite counterterm 

- i J dVrA (B.24) 



15 In making the comparison one should keep in mind that R\jkl" = — nyjy 
two typos in (3.12) of [33]. Namely, in the third line of (3.12), — g(2)g(2)ij should read Tr g (2)5(2)ij) while in the 
beginning of the fourth line RikRj should read 2RikRj. 
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